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I. 



INTRODDCTION 



A. PROBLEM STATEMENT AND OBJECriVE 

The prediction of meridional flows within turbo machines, 
be they compressors or turbines, is a difficult but 
important part of the design process. The difficulty arises 
from the presence of three- dimensional and viscous effects 
within all turbomachines and the importance arises from the 
necessity to design accurately and efficiently. 

To simplify the problem of viscous, three-dimensional 
analysis, Wu [ Ref. 1 ] showed that this complicated flow may 
be analyzed by solving two interrelated flows: one being the 
blade-to-blade flow describing the flow between rotating 
blades and the other being the meridional through flow which 
describes the radial equilibrium. These flows are depicted 
in Fig 1. In addition, an inviscid and axi symmetric 
assumption is made in the through-flow thereby simplifying 
the flow to a two- dimensional, axi symmetric, inviscid, and 
compressible analysis. 

Three methods may be found in current reports regarding 
the solution of the radial equilibrium equation. The first 
two are the streamline curvature method [Ref. 2, 3, and 4] and 
the matrix method [Ref. 5 and 6] which is oasically a finite 
difference technique. The third, a relatively new method, is 
the finite element method. As shown by Hirsch and Narzee 
[Ref.7]', the solution of the radial equilibrium equation by 
the finite element method is achieved by arranging the 
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the radial equilibrium equation is presented followed by the 
application of the finite element method to this equation. 
Section III describes the computer program in some detail. 
Section IV contains selected test cases which were used in 
program testing and checking. In Section V, conclusions are 
presented along with recommendations for further study and 
work on the project. The appendices conrain the program 
listing along with a sample test case for reference by the 
user. In addition, a list of references is contained for 
further reading on the subject of this paper. 
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I:E: 




Figure 1 - MERIDIONAL AND BLADE- TO- BLADE PLANES 
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II. 



THEORY 



I 

|a. the derivation of the radial squilibrioh equation 

The following discussion is taken from Ref. 7 with 
slight changes in notation. The basic turbomachine geometry 
^to be analyzed is depicted in Fig 2. Although the machine 
noted is one stage of a compressor, a similar analysis to 
|the one that follows may be applied to other machines such 
as axial turbines and mixed-flow machines. 

One begins with the Euler equation assuming the viscous 
forces to be negligible. 

^ + (v- v) v’ - vp|f (II. A. 1) 

The continuity equation, assuming unsteady flow is, 

^ + v(fv’)=0 (IX. A. 2) 



The First Law of thermodynamics in a fluid field 
jecomes. 



Tvs = vh - vplf 
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OUTER CASING 




AXIS OF SYMMETRY 



Figure 2 - TURBOMACHINE GEOMETRY 
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\ 



substituting equation (II. A. 3) 
leads to the Crocco equation. 



into equation (II. A. 1) 



^-v4vxv)= Tvs - VH 



(II. A. 4) 



where H is the total enthalpy- 

■ , a steady ana adiabatic flow, the energy 

Assuming u snea y 



eq 



uation becomes simply. 



(v-V^H=0 



(II. A. 5) 



which shows that along a streamline 
the total enthalpy is constant. 



iu a stationary 



system. 



In a relative system, such a^ 
row, the total relative velocity, W, 
following form. 



the case 
can be 



in a rotor blade 
expressed in the 



—■7 

V/ 



- y ou> A & ” 



—7 



^■7 -■> 

V + u. 



(II. A. 6) 



.hare « is the constant angular velocity and d is the 
constant peripheral speed of the relative syste*. 

»cv. the Crocco eguation in a relative systee becones, 

— I (II. A. 7) 

<^"t 

0 .- ITT 1 ’=,) for the stationary system, the 
parallel to equation (II. A. 5) f /relative) 

• and adiabatic ^rexa i 

energy equation, assuming steady 

flow in a relative system, becomes 
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$ 



(w’- 7 )Hj =o 



(II. A. 8) 



where H is the relative total enthalpy expresed as 




From the following velocity diagram, 




follows , 
(II. A. 9) 
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equation (II. A. 9) may be arranged as follows. 



Since, 



w 







-- -vUs -U)" 



(II. A. 10) 



then , 



VJ% - ZUYe +U"' 



(II. A. 1 1) 



and. 



W^“ v"*'+ ~ 



(II. A. 12) 



Substituting equation (II. A. 12) into equation (II. A. 9) 
leads to the following relation. 



h * X ' Ms = H-Ws (II. A. 13) 



Equation (II. A. 8) shows that is constant along a 

streamline in a relative system. 



Upon integrating equation (II. A. 8) between the rotor 
inlet and outlet, the Euler equation for turbomachines is 
found , 



AH 



OU-i 



ifj 



A (U'V&j 



(II .A . 1 4) 
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It may be shown [Ref. 9] that by circumf eren tially 
averaging equation (II.A.1), and under the axi symmetric 
flow assumption the following relation is valid, 

-Vx(v^^)-T-VS-VH (II. A. 15) 

where is the body force of the blades acting on the fluid 
and all variables are mean values along the direction of the 
circumference. Hence equation (II. A. 15) is an approximation 
for axi symmetric flow. As a final note on equation 
(II. A. 15), since the viscous forces were neglected in 
equation (II. A. 1) , there must be a force introducing the 
entropy variations along the blade. This force is 

proportional to the pressure loss coefficient and is labeled 
Fj , the dissipative force. F^ produces work which in turn 
produces entropy production radially along the blade. Under 
the axi symmetric assumption, entropy varies axially and 
radially only and is assumed to be proportional to the 
pressure loss coefficients [Ref. 2 and 8]. 



* Due to boundary conditions imposed on the problem and 
the axi symmetric assumption, cylindrical 
coordina tes , (r ,0, z j , will be used in all subsequent analysis. 
Therefore, equation (II. A. 15), in cylindrical coordinates 
and with axial symmetry is as follows. 



^ ^ (p\j ] vl^Vz-A-v] — 












(II . A. 17) 



(ir . A. 18) 






iH rAi 

it ■ ' 




It is important to note here that under the axi 
symmetric assumptions, eguation (II. A. 15) reduces to the 
following. 



V' 




o 



(II. A. 19) 



Likewise in a relative system (rotor) , the axi symmetric 
assumption leads to the following. 



W- F\ ^ ^ (II. A. 20) 



Equation (II. A. 16) describes the meridional through flow 
radial equilibrium equation for the finite element method. 
Since one is concerned with the meridional plane, the 
following derivative expression is taken from Fig 3. 






(II. A. 21) 



Therefore equation (II. A. 17) reduces to. 






(II .A . 22) 



which reveals that in a duct where there are no blades and 
therefore no blade forces, angular momentum is constant 
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along a streamline. In that case, 




(II. A. 23) 



As shown in Ref. 9, the circumferentially averaged 
continuity equation is the following, 

(II. A. 24) 

where b is the blockage factor defined by Hirsch and Warzee 
as the tangential area reduction due to the thickness of the 
blade. 




(II. A. 25) 



where t is blade thickness and s is blade spacing. 
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Figure 3 - MERIDIONAL PLANE 
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One further step in the formulation of the radial 
equilibrium equation for solution by the finite element 
method involves introducing the stream function. In 
cylindrical coordinates, the stream functions are defined as 
follows , 




I ^ 

'I 



(II .A .27) 



Substituting these expressions into equation (II. A. 16), 
the equation becomes. 



Jl 





(II . A. 28) 



The right hand side of equation (II. A. 28) is applicable 
to the absolute flows in the stator and duct regions. For 
relative flows such as those in the rotor, the right hand 
side is modified by replacing the total enthalpy, H, by the 
relative total enthalpy, H, and the quantity 7e/R is 
replaced by We/R. 



As a last assumption in the formulation of the governing 
relation for the meridional through flow radial equilibrium 
equation, both the radial component of the body force, 
and the radial component of the dissipative force, F^ , are 
neglected. This assumption, [Ref. 1,8] does not hamper the 
accuracy of the results for conditions at design speed. 
Even though published compressor performance data used for 
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the test case in this thesis was obtained at 0.5 design 
speed, these force terms were also neglected in the computer 
program. As will be shown later, this assumption could 
possibly have had adverse effects on the predicted axial 
velocity profiles at the rotor hub and tip regions. 



The final representation of the meridional radial 
equilibrium equation to be solved by the finite element 
method is as follows. 



where , 



A f, A ^ 

Vk ^ J >% 









= o 



(II .A. 29) 



(II .A. 30) 



and 



t - — 1 (pm ) 



(II. A. 31) 



B. THE FINITE ELEMENT METHOD APPLIED TO THE RADIAL 
EQUILIBRIUM EQUATION 
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(II. B. 1) 



i. 




O 



where this equation represents the flow in the volume, V. 



The boundary condition for this partial differential 
equation, after dividing through by R, is, 

CH'-i)! =o 

where this equation solves the flow on the closed boundary 
of the volume, or, S. 



By applying the weighted residual process to equations 
(II.B.1) and (II. B. 2) and using an arbitrary weighting 
function, W(r,z), one has 







(II. B. 3) 



where r and r, , are the volume and surface residuals 

VOL 

respectively, or. 





Q (II. B. 4) 






(II. B. 5) 



If the solution to equation (II.B.1) was exact, both r^^^ 
and would be equal to zero. 



In order to clarify the boundary condition, equation 
(II. B. 2) , one may analyze the equation as follows. 
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On the surface ,S , where ^ 



is specified. 




(II. B. 6) 



and. 



<=o 



(II. B. 7) 



Similarly, 



on the surface. 



where ^y\ ' ^ , 



where 



(II. B. 3) 

^ 1^0 



S'.n =o, US', ^ S' 



(II. B. 9) 



Due to the axi symmetric assumption, the final equation 
will not involve dV and dS but the intersection of dV and dS 
with the meridional plane. Therefore, one must transform the 
volume integral, dV, to a surface integral and the surface 
integral, dS, to a line integral. 

Hence, let, 

djl = intersection of dV and meridional plane 
dC = intersection of dS and meridional plane 

J .S' ^ 2TT e d C. 

With this transformation, one may rewrite equation 
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(II. B. 3) as follows. 




where on the contour, 



(II .B. 10) 



One must now integrate the first term in 
(II.B.10) by parts to obtain the following, 

- [w-j-cl-rz, 

-n. 



eg nation 



(II .B. 11) 




Inspecting the first term in equation (II.B.11) , one may 
use the following integral theorem to simplify further 




(II. 3. 12) 



Rewriting equation (II.B.11) , gives, 

-h 2 t7c1C O 

c 

Finally, since 



(II .3. 1 3) 








(II .B. 14) 



equation (II.B.13) reduces to the following. 
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(II. B. 15) 



One now has the final equation in the form for use by 
the weighted residual method using any arbitrary weighting 
function, W (r,z) . As noted previously, the Galerlcin's 
Method will be used here which implies that the weighting 
functions are the same functions used in approximating the 
stream function, ^ . 

Before applying the finite element method, one must 
discretize the continuum and then approximate the unknown 
function , by a set of polynomials. For this particular 
problem, eight- noded iso-parametric elements were chosen for 
discretization, see Fig 4, and the following approximating 
functions were used. 



location within the element. The shape functions, , used 
here are defined by the following relations as shown in 
Ref . 10, 




(II .3. 16) 



where , 



~ shape functions 
= value o f 'f’ at the node 
= value of at a 



at any arbitrary 
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(II .3. 1 7) 



Mu)- 

i(i* ('-'(’■) 

where the following coordinate transformations are used 

S 







I 



. 1 8 ) 
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Figure 4 
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V|^ = +1 




5 6 7 






= -f-i 
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= -1 



- ISOPARAMETRIC QUADRILATERAL ELEMENT 
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At this point, one is ready to apply the Galerkin Method 
to equation (II.B. 15) by substituting equation (II. 3. 16) for 
the unknown function*^ , and N- for the weight function, W, 
which yields. 




Q (II . B. 19) 



This integration yields the following system of 
equations which is solved for the unknown nodal'Vf 



where. 




(II.B. 20) 



and , 





XH. XU: 
>r 






(II.B. 21) 



-A 



(II.B. 22) 



In addition since both the 'stiffness matrix', K, and the 
right hand side vector, [F], are functions of , the system 
as defined by equations (II. 3. 20) through (II. 3. 21) must be 
solved iteratively. 

At this point one has the total finite element 
formulation of the radial equilibrium equation as defined by 
equations (II.B. 19) and (II. 3. 20). The problems which 
remain to be clarified are basically two fold. Firstly one 
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must evaluate the integrals in equations (II.B.19) and 
(II. B. 2) by numerical methods, and secondly, the solution 
procedure for the non-linearity must be formulated. In Part 
C, both of these final steps are presented. 

C. NUMERICAL INTEGRATION OF STIFFNESS MATRIX AND SOLUTION 
PROCEDUR E 



1 . Nume ri c al inL§3E5.ii2S 21 2]l2 stiffn ess matrix 



(II.B.21) must be perfacmed numerically. In addition, one 
realizes that the derivative expressions enclosed within the 
interval must be evaluated by a coordinate transformation. 
This is done in the following way. 



As noted in Section II. B, evaluation of equation 



Since 
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(II. C. 1) 



S 






then 



■ 







(II. C. 2) 
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and in matrix form. 





(■ 


H ir'' 






. = 1 










W ir 








in nl 





(II. C. 3) 



Furthermore, defining the Jacobian matrix as. 




(II. C. 4) 



then by dividing both sides of equation (II. C. 3) by J, one 
has the following transformation. 




(II. C. 5) 



In addition, it has been shown [Ref. 9] that 



JiAr UU^dv| 



(II. C. 6) 



Now, with equations (II. C. 5) and (II. C. 6), equation 

(II.B.21) becomes the following, 
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(II. C. 7) 



Equation (II. C. 7) is best integrated using the 
Gauss-Legendre integration method since it is of the 
following form. 



' I 





(II. C. 9) 



where A. and B^are coefficients (Fig 5) for both two and 
three point Gaussian Quadrature. 



At this point, one has the tools to calculate all 
the element s of the stiffness matrix. In like manner, the 
right hand side vector, f, is calcualted by numerical 
integration . 
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Figure 5 - GAUSSIAN INTEGRATION POINTS 
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• Soluti on pro ce dure 



The following is a synopsis of the basic solution 
process. Specific details concerning equations and methods 
of computer coding are covered in the proceeding section. 
The proceeding is meant to give the reader a preview of the 
solution process. 

a. Discretization 



Initially the machine un 
discretized into eight-node iso-parame 
axial calculation stations are placed a 



duct regions 


and along 


blade edges 


and ce 


and stator as 


shown in 


Fig 6. At 


this 


topology and 


nodal coor 


dinates are 


specif 



der analysis is 
trie elements. The 
rbitrarily in the 
nters for the rotor 
point the system 
ied . 



b. Initialization 



To begin the iteration process, one m 
an initial internal stream function, velocity, an 
distribution. In the program, the initial inter 
function was assumed to be that of the outer 
throughout while the velocity and density distri 
assumed to be that of the inlet. 



ust 
d d 
nal 
bo 
buti 



assume 
ensit y 
stream 
undary 
on was 
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Figure 6 



COMPRESSOR 



DISCRETIZATION 



33 



c. Calculation of thermodynaraic variables 

Before calculating the right-hand side vector, f 
, one must obtain distributions of angular momentum, 
enthalpy, and entropy. This is done by first calculating 
the thermodynamic variables at the inlet axial station from 
the given inlet conditions. In order to proceed axially 
through the machine to calculate the nodal angular momentum, 
enthalpy, and entropy, the following three equations derived 
in Section III are used. 



H ^ C|>T 



constant along a stator streamline 




constant along a rotor streamline 
constant along a duct streamline 



An example of this calculation procedure for the duct region 
is shown graphically in Fig 7. 
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Figure 7 - DUCT ELEMENT 
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In this figure, the angular momentum at point X 
is equal to the angular momentum at point Y. Wore formally. 




(II. C. 2. 1) 



Since the previous axial station's thermodynamic 
‘variables are known, one must now find the values of ^ and 
at point Y. This is done iteratively in the following 
way. Since 









(II. C. 2. 2) 



and along the left side of the element. 




(II. C. 2. 3) 



then equation (II.C.2.2) may be solved for by a suitable 

iteration method. As will be shown in the next section, a 
half-interval method was used to obtain the unknown >| 
Once is known, then equation (II. C. 2.1) is solved for 

the angular momentum at point X. The rotor and stator are 
handled in a similar fashion. In addition, the rotor and 

stator deviate the flow creating a three-dimensional flow 
field between the blades in the respective blade row. Low 

speed cascade correlation data [Ref. 13] was used to 

calculate the effective turning angles in the rotor and 

stator. These effects are calculated beforehand with known 
mass flow rate and uniform axial velocity assumptions at the 
rotor inlet. The results of these calculations are part of 
the input data routine in the form of relative and absolute 
flow angles at the rotor nodes and absolute flow angles at 
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the stator nodes. This will be shown more exactly in the 
next section. 



d. Calculate matrices 



At this point, the right hand side vector, f, and 
the stiffness matrix ,K, are calculated. 



e. Solve system of equations 



The system of equations as shown in equation 
(II. B .20) is solved for the nodal stream function. 

f. Perform relaxation iteration 

Due to the strong non-linear properties of the 
system of of equations, the following iterative scheme is 
necessar y. 

4 oC (II.C.2.4) 

t I i, _ 

where c( is the under relaxation factor. As will be shown in 
Section III, this scheme is performed only in certain 
regions of the machine and in addition after a specified 
number of iterations. 

g. Update velocity and density profiles 

Using the current nodal distribution of the 
stream function, axial and radial nodal velocity components 
are calculated along with a new nodal density distribution. 
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¥ 

I 



Again, this calculation procedure will be shown in the next 
section . 



h. Test for convergence of ^ 



Stream function convergence criteria is now 
tested and will determine if further iterations are 
necessary. The solution is said to converge if the following 
equation holds for all nodes. 

, n +- 1 



't;- 'v; 






Y1 + ' 






(II. C. 2. 5) 



where € is a designated requirement for convergence. 



i. Summary 



In summary, the eight steps involved in the 
solution are noted below; 

(1) Discretize the continuum. 

^ 2) Assume an initial stream function, velocity, 
and density solution. 

(3) Calculate the nodal thermodynamic variables 
from the given inlet conditions. 

(4) Form the right hand side vector, f(r,z), and 
the stiffness matrix, K. 

(5) Solve the system of equations, given by, [K] 
= [F] for a new stream function distribution. 
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(6) Perform relaxation iteration if required, 

(7) Calculate new nodal velocity and density 
distributions from the current stream function solution. 

(8) Test the solution for convergence, and if 
required, repeat steps (3) through (8) using the current 
nodal stream function values. 

This concludes the solution description and now 
one is ready to more completely understand the computer 
program which assembles the preceding eight steps. 
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III. THE PROGRAM 



A. OVERALL FLOWCHART AND DESCRIPTION 



The overall flowchart of the program is depicted in Fig 
8. Those blocks denoted by the letter ’S’ are subroutines, 
while the remaining calculations are an integral part of the 
main program. 

After proper dimensioning of all arrays and subsequent 
initialization, the input data are read and then printed. 
This not only presents a physical picture of the problem but 
also serves as a cross check to the user for correct data 
insertion. In addition, a subroutine is available to obtain 
a computer drawn plot of the mesh (Fig 6) and is a further 
check on proper data input. 

At this point all the necessary variables have been 
stored and the iteration counter for stream function 
convergence is set. With the current nodal values of and 
the given inlet thermodynamic conditions, the thermodynamic 
variables throughout the machine are calculated. From the 
calculated values of enthalpy and angular momentum, 
(isentropic flow is assumed) , the right-hand side vector is 
calculated followed by the stiffness matrix calculation 
(equation II . B .2 1 ) . 

The system of equations (equation (II.B.20)) is now 
solved for the new nodal stream function distribution. It 
is here where for all iterations but the first that a 
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relaxation factor is applied as noted previously in equation 
(II.C.2.4). The reasoning behind not applying the 
relaxation scheme to the value of nodal after the first 
iteration is the face that the first iteration produced a 
close approximation to the correct stream function 
distribution. With this close approximation to the stream 
function came a velocity and density distribution which in 
turn was near the correct solution. It was found that if 
the first iteration was relaxed, the second iteration became 
unstable since in fact the velocities and densities were 
themselves farther from the true values than were assumed 
initia 11 y . 

After testing the nodal stream function for convergence 
by use of equation (II.C.2.5), the calculation process is 
either repeated or ceased by virtue of convergence or 
limiting the number of iterations. 

As stated previously, low speed cascade correlation data 
[Ref. 13] were used to calculate turning angles in the blade 
regions. These angles were assumed constant throughout the 
solution and not refined after subsequent iterations. 
Further work on the computer prog ram could entail an 
additional computational routine which would calculate the 
new turning angles after each iteration. A sample 
calculation of rotor turning angles is shown in Appendix D. 

In the following sections the program structure is 
examined in more detail. 



4 1 
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B. THE MAIN PROGRAM 



1 . input routine 

The following is a description of the inp 
required by the program. The data are arranged into 
categories described in the following manner. 

a. category 1 

Problem identification. 



b. category 2 



Number of nodes and number of elements. 



c. category 3 

Node numbers, nodal coordinates and 
blockage factor. 

d. category 4 
System topology. 

e. category 5 



t data 
tvel ve 



nodal 
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Element type; duct, rotor, or stator. 



f. category 6 

Absolute flow angles for rotor and stator nodes. 

g. category 7 

Relative flow angles for rotor nodes. 

h. category 8 

Inlet thermodynamic quantities. 

i. category 9 

Physical constants for fluid under observation. 

j. category 10 

First estimate of internal stream function. 

k. category 11 

Node numbers and specified nodal stream 

function . 

l. category 12 
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Node numbers where the right hand side, f(r,z), 
is to be calculated. 

Before describing in detail the format to be 
followed for data insertion, it is important to note the 
following assumptions. 



outlet . 



(1) Uniform flow conditions at inlet and 



(2) Uniform flow conditions at rotor inlet 
for calculation of appropriate turning angles. This 

assumption is necessary to calculate the values of rotor and 
stator flow angles. 

With this in mind , the discussion will 

CO n t i n u e . 

The following describes each category in mors 

deta il. 



Category 1; 

Format: (20 A4) 

Number of cards: 1 

Procedure: Enter the title of the problem in 

columns 1-20. 

Category 2: 

Fo rm at ; (2110) 

Number of cards: Equal to the number of 

nodes in the system. 
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Procedure: Enter the number of nodes in 

columns 1-10, and the number of elements in 11-20. Both 
integers must be right justified. 

Category 3: 

Format: (I10,3F10.0) 

Number of cards: Equal to the number of 

nodes in the system. 

Procedure: Each card contains the node 

number followed by the 2 coordinate, R coordinate, and nodal 
blockage factor. The coordinates are in dimensions of 
inches . 



Category 4: 

Format: (9l5) 

Number of cards: Equal to the number of 

elements in the system. 

Procedure: Each card contains nine integers 

right justified in columns 5, 10, 15, etc., through 45. The 

first integer is the element number followed by the eight 
nodes associated with that element. It is important to note 
that the nodes are read in starting with the upper right 
hand node and proceeding in a counterclockwise fashion 
around the element. 

Category 5: 

Format: (2110) 

Number of cards: Equal to the number of 
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elements. 



Procedure; Enter the element number in 

columns 1-10, followed by the integer 'V (duct), '2* 

(rotor), or '3* (stator) describing the element as either in 
a duct, rotor, or stator region. 

Category 6; 

Format: (6X, A4 ,I 10,F10 .0) 

Number of cards: Equal to the number of 

rotor and stator nodes plus one 'STOP' card. 

Procedure; Enter the node number (right 

justified) in columns 11-20 followed by the value of the 

associated absolute flow angle in radians in columns 21-30. 
The last card in this category is a 'STOP' card entered in 
columns 7-10. 



Category 7: 

Format: (6 X , A4 ,I 1 0 , F 1 0 . 0) 

Number of cards: Equal to the number of 

rotor nodes plus one 'STOP' card. 

Procedure: Enter the node number (right 

justified) in columns 11-20 followed by the value of the 
associated relative flow angle in radians in columns 21-30. 
The last card in this category is a 'STOP' card. 

Category 8; 

Format: (7F10.0) , (F10.0) 
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Number of cards: 2 



Procedure; Enter the following quantities in 
the prescribed order and with the noted dimensions. 

First card 

Mass flow rate: (Ibm/sec) 

Inlet axial velocity: (ft/sec) 

Outlet axial velocity: (ft/sec) 

Inlet total density: (Ibm/ft^ ) 

Inlet static density; (Ibra/ft^ ) 

Inlet total pressure: (Ibf/in’’ ) 

Inlet total temperature: (*R) 

Second card 
Speed: (RPM) 

Category 9: 

Format: (3 F 1 0 . 0) 

Number of cards: 1 

Procedure: Enter the following quantities in 

the prescribed order. 

Gas constant: ( f t- Ibf /Ibm- *R) 
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Ratio of specific heats 



Co nstant 

(BT0/lbm-*’R) 



pressure specific heat: 



Category 10: 

Format.: (FI 0 . 0) 

Number of cards: 1 

Procedure: Enter the first estimate of the 

internal stream functiion to be used in the first iteration. 

Category 11: 

Format: (5X,A4,I10,F10.0) 

Number of cards: Equal to the number of 

nodes having a specified value of the stream function plus a 
•STOP' card. 



Procedure: This set of cards allows the 
stream function boundary conditions to be read in. A 
typical card contains an integer, right jusrified in columns 
11-20 , which is the node number, followed by the value of 

the specified stream function in columns 21-30. The last 
card is a ’STOP’ card. 

Category 12: 

Format: (5X, A4 ,1 10) 

Number of cards: Equal to the number of 

nodes where the right hand side is to be specified. 
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Procedure: Enter the node 
justified in columns 11-20, where the right 
be calculated. Again, the last card in this 

’STOP' card. 



number, 
hand side 
category 



right 
is to 
is a 



After all the data 
the input data is printed 
verification by the user. The 
Appendix C. 



has been read by the program, 
and the mesh is plotted for 
sample format is shown in 



section 

K. 



This concludes the input 
describes the calculation of 



routine. The next 
the stiffness matrix. 



2 . Stif ^ess matrix evaluation 



As shown Dceviously in section II.C.1, the following 
eguatlon describes each ters in the eight by eight eleeental 

matrix. 






(III. 3. 2. 1) 



in addition, 'h' is defined in the following way in order to 
numerically integrate the equation. 

1 

K ^ — I (III.3.2.2) 






where b is defined as the elemental blockage factor taten as 

an awerage over the eight nodes of the particular element 

and (VV are the defined 3 auss-,uadrature integratron 
points . 
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Figure 9 - STIFFNESS MATRIX EVALUATION 
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Fig 9 depicts the flowchart for both elemental 
stiffness matrix evaluation and the assemblage into the 
system stiffness matrix. More specifically, the figure 
shows a three-point Gaussian Quadrature scheme but can be 
changed to a two-point scheme by simply integrating four 
times instead of nine as shown. 



The actual coding of the stiffness matrix evaluation 
and assemblage may be found in lines STR03510 through 
STR04770 in the computer program. 



3 . Solution of systems of §auat ions 



At this point, the system of eguati 
for the boundary conditions and solved for 
function values. An equation solving 
available in the system library was used 
It was found that no comparable savings 
using a banded equation solver. 



ons are modified 
the nodal stream 
routine, DSIMQ, 
for this purpose, 
was realised by 



4, Iteration schemes 



As noted previously in Section II. C, a relaxation 
scheme is necessary for convergence to a solution. 

Two distinct differences with regard to the 
iteration method were noted from that of Ref. 7. Firstly, it 
was found that relaxation was necessary only in the rotor 
and stator elements and also in the duct region between the 
rotor outlet and stator inlet. Secondly, due to the extreme 
non linearity in the rotor-stator areas, a switch was 
required which changed the sign of in equation (II. C. 2.4) 
as required for stability of convergence. Claritication of 
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this change follows: It was found that during the initial 
three or four iterations, the stream function values of the 
rotor-stator nodes sometimes exceeded the value of the upper 
boundary. Due to an abscence of sources within the domain 
of solution, this occurrence was incompatible with the 
boundary conditions. At this point, it was necessary to 



make d. 


negativ 


e in equation 


(II. C. 2.4). 


During subsequent 


iterations, as 


the solution 


converged. 


the rotor-stator 


regions 


became 


stable and the 


sign of d. 


was returned to 


its 


positive 


val ue. 


This iteration proved 


to stabilize 


the 


solution 


w ith 


respect to 


stream fu 


notion values 


and 


velocities. 











The iteration procedure is coded in the computer 
program from lines STR05070 through STR05200. 
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C. THE SUBROUTINES 

The following describes each of the six subroutines in 
the computer program. Each subsection contains a list of 
calling arguements and for subroutines ?CAL, SEINE, and VEL, 
a basic flowchart. In addition, for those subroutines whose 
mathematical theory was nor presented in Section III, a 
brief treatment is also given. 
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^ • Su broati ne sh ape 



This subroutine calculates the shape functions 
(equation (II.B.17)) at the values of ^ and as requested 
in the arguement list below. 

SUBROUTINE SHAPE (E,Z,SF) 

E = value of ^ (input! 

Z = value of (input) 

SF = eight by one vector of the eight shape functions. 



2 . Subroutine qacob 

JACOB calculates the Jacobian matrix as defined in 
equat ion (II . C. 1 . 4) for the value of ^ denoted in the 

arquement list. 



SUBROUTINE JACOB ( E 1 , Z 1 , D , S , RC $, Z C$ , RJ AC) 

El = value of (input) 

Z1 = value of ^ (input) 

D = eight by one vector of (calculated) 

E = eight by one vector of (calculated) 

RC$ = eight by one vector of the 'r' coordinates of the 
nodes associated with the element (input) 

ZC$ = eight by one vector of the 'z' coordinates of the 
nodes associated with the element (input) 

RJAC = two by two Jacobian matrix (output) 
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In addition, the subroutine assumes that the vectors 
RC$ and ZC$ contain element coordinates arranged in a 
counter clockwise fashion beginning with the upper right 
corner node. 



3 • S ubroutine si in e 

This subroutine calculates the thermodynamic 
variables throughout the machine given the inlet conditions 
as described in Section II. C. 2. The calling arguements are 
defined below. 

SUBRODTINE SLINE (BIN LET, RC , P SI , WR L , H , D7EL, V VE L ,T VE L, 
NODE, NNO DEI, CP,TT, KK , ALP , tf G , TWEL, BE , HS) 

UINLET = Inlet axial velocity 

RC = Nodal ' r' coordinates vector 

PSI = Nodal stream function vector 

WRL = Nodal angular momentum vector 

H = Nodal total enthalpy vector 

UVEL = Nodal axial velocity vector 

VVEL = Nodal radial velocity vector 

TVEL = Nodal absolute tangential velocity vector 

NODE = Matrix containing nodes associated with the 
element 

INLET = Vector containing node nu-mbers at inlet station 
NNODSI = Number of nodes at inlet station 
CP = Specific heat 
TT = Total temperature at inlet 
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KK = Iteration counter 

NTE = Element type vector 

ALP = Nodal absolute flow angle vector 

TWEL = Nodal relative tangential velocity vector 

BE = Nodal relative flow angle vector 

HS = Nodal static enthalpy vector 

As shown in Fig 10, the basic calculation procedure 
begins with calculating the required energy and momentum 
values at the inlet station. At this point, beginning with 
element one, the element type is interrogated to distinguish 
between duct, rotor, and stator elements. If the element is 
in a duct region, then the streamline intersections for 
local nodes 2,6, 7, 3 and 1 (Fig 7) are determined along with 
the associated values of energy and angular momentum. For 
the rotor and stator elements, one must initially find the 
energy and momentum values at local nodes 3,4,5 (Fig 7) due 
to the discontinuities imposed by the blade edges. Once 
these calculations are performed, then the process for the 
remaining nodes in the element proceeds in a similar fashion 
to the duct elements. 



After all the 
new distributions of 
returned to the main 
Specifically, these 
subroutine , FCAL, for 
vector . 



elements have been cycled through, the 
nodal angular momentum and energy are 
program for further computations, 
values will be used by the next 
calculation of the right hand side 
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Figure 10 - SUBROUTINE SEINE 
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4 



• Sub routine fcal 



FCAL 

defined by 
ident ica 1 
integration 
to be coded 

1 

T ^ 



calculates the right hand side vector as 
equations (II. A. 31) and (II.B.21). Using the 
coordinate transformations for numerical 
as described in Section II. C, the final equation 
is the following, 

ffJL 

Ihm, < 

(III. C. 4. 1) 






where isentropic flow is assumed, and. 



Wj^ = angular momentum 

(III.C. 4. 2) 



H; = total enthalpy 

V 



The arguement list is defined below. in addition, 
only those variables in the list, which have not been defined 
previously are described. 






SUBROUTINE FCAL (F,M,H,ZA,SA,UVEL,RC,ZC,NRL 



NODE, NN 
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TW 
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is to be specified 

NN = number of nodes 
NE = number of elements 

NNFSP = Number of nodes where the right hand side is 
specif ied 



Fig 11 depicts the basic flowchart for the 
subroutine. To initialize the procedure, one begins with 
the first node (upper right hand corner) of the first 
element. A switch is then applied which determines if the 
right hand side is to be calculated at the node or if a 
stream function value has been specified. This information 
is transferred from the main program through the arguement 
list. Once the node is allowed through the switch, then the 
integration process is started at the first integration 
point. As in Section III. A. 2, the flowchart depicts a 
three-point Gauss Quadrature scheme. After the integration 
has been completed, a switch determines if all the local 
nodes in the element have been cycled through and if so, 
then the assembly of the elemental vector, F$, is performed 
to build the system right hand side vector , F. Finally, 
the subroutine^determines if all the elements have been 
examined in order to signal completion of the right hand 
side vector. At this point, the vector, F(r,z), is returned 
to the main program for problem solution. 
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Figure 11 - SUBROUTINE FCAL 
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5 . Subr out ine vel 



This subroutine calculates axial and radial 
velocities and also densities at each of the nodes from a 
known stream function distribution. As noted previously in 
Section II. C. 2, both velocity and density profiles are 
updated after obtaining the latest value of nodal stream 
function . 

The velocity calculation proceeds from the stream 
function equations. 









I ^ 



(III. C. 5. 1) 



(III. C. 5. 2) 



where 'b' is the tangential blockage factor. Since r, ,and 
are of the following form, 

t 






1- i 

L •' 



(III. C. 5.3) 



then the equation for the axial velocity, , becomes. 



V,- 



i 















(III.C. 5. 4) 



Again, since the shape function, (^/''|^) / is not 



an 
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implicit function of 'r' and 'z', one must use equation 
(II.C.1.5) to obtain the proper derivatives for computation 
of equation (III. C. 5. 4). For example, from equation 
(III.C. 5) , 




(III . C. 5. 5) 



At this point, with equation (III.C. 5. 5) 
substituted into equation (III.C. 5. 4) , one has the complete 
expression for the axial velocity as functions of^ One 
proceeds similarly for expressing the radial velocity, 7 , 
in terms of ^ andY^. 

In order to calculate the nodal density, one uses 
the following density relation for flows in the stator and 
duct regions. 

1 - ('i- IdwA'M 

^ ^ (III.C. 5. 6) 



where 






is the stagnation density. 



Since , 



(Vj . + tM.V) ( I II . C. 5 . 7) 



then , 



1 . 




tU /__L 




J_ 



(III.C. 5.8) 
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since the density appears on both sides of the 
equation, the new nodal density is obtained iteratively at 
the node. 



For the relative flows in the rotor, the following 
relation for static density is used [Ref. 14]. 



1 




2 - d.' 



r- 



(III.C. 5.9) 



Again, the solution of the nodal density is obtained 
in an iterative fashion. 

In the following arguement list, only those 
variables not defined in the previous subroutine 
descriptions are noted. 



SOBRCOTINE VE L ( NE , NN , RC , NODE , G , RG , TT , RKOT , RHCN , ZC, 

PSI,RHO, B,UINLET,DVEL,VVEL,RHOSTA, NTE, ALP) 

G = Ratio of specific heats 

RG = Gas constant 

RHOT = Total density at the inlet 

RHON = Wor)c vector which contains the new nodal density 
distribution 

RHO = Nodal static density vector 

B = nodal blockage factor vector 

RHOSTA = Static density at the inlet station 



The basic flowchart for SUBROUTINE VEL is shown in 
Fig 12. Beginning with the first node of the first element, 
the Jacobian matrix (equation (II.C.1.4)) and its inverse 
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I 



are found. At this point the partial derivatives with 
respect to *r' and ’z* of the shape functions are found as 
noted in equation (III. C. 5. 5). A switch then allows those 
nodes not at the inlet station to pass and calculates the 
new density and velocities at the nodes. For those nodes at 
the inlet, the velocities and static densities are retained 
at the given inlet conditions. This is done to maintain 
boundary condition integrity for the solution. After 

cycling through all elements, the subroutine returns the new 
nodal velocity and density distributions to the main 
program. 
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Figure 12 - SUBROUTINE VEL 
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6 show the computer predictions plotted with the 
values and the numerical solutions obtained by 
arzee. The profiles shown were obtained after 
ons and using a relaxation factor of 0.2. The 
w that the best overall agreement with 
data occurred in the stator inlet and outlet, 
on the worst error was 17% which occurred at the 
inlet. The average error throughout the stator 
respect to experimental data was 6.6%. 
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The rotor hub and tip outlet area exhibited 
instabilities in density convergence using equation 
III. C. 5. 9. Specifically, the density solution converged to 
within 8^ at the rotor outlet tip and hub. It was found 
that by not allowing the nodal density at these nodes to go 
below a critical value of 0.06 Ibm/cu ft, the solution for 
the stream function converged. By allowing the nodal 
densities at the rotor outlet tip and hub to go below this 
critical value, the computed velocities at these nodes 
became increasingly large and the arguement within the 
brackets of equation III. C. 5.9 became less than one. This 
prevented continuation of the iterations for the stream 
function solution. In addition, the rotor tip outlet 
exhibited more instability than the rotor hub outlet. The 
static density at the rotor hub outlet oscillated about a 
value of 0.062 Ibm/cu ft while the rotor tip outlet was 
constantly driven to the critical value of 0.06 Ibm/cu ft. 
One method attempted to alleviate this problem was the 
following. Since a hal f- interval iteration routine was 
used, one trial run involved reversing the direction of 
consecutive guesses when the density iteration did not 
converge. It was found however, that after three to four 
iterations of the system of equations, the static densities 
at the rotor outlet tip and hub were again driven to smaller 
and smaller values which led to instability once more. The 
nodal densities converged at all interior points of the 
rotor edge and mid-blade regions and also at all the rotor 
inlet nodes. By including all rotor nodes, the average 
error with respect to experimental data was 27.5%. 

Fig 17 shows a plot of convergence criteria, £, versus 
the number of iterations for a relaxation factor of 0.2. 
The stability of convergence is shown to initially decrease 
and then after the third iteration oscillates about an 
approximate value of 28%. It is important to note that this 
curve represents the maximum value of 6 as shown in equation 
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CONCLUSIONS AND EICOMMDillQNS FOR FURTHER STUDY 
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Figure 14 - AXIAL PROFILE AT ROTOR OUTLET 
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APPENDIX A 
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NASA TASK-1 TRANSONIC COMPRESSOR 
1C7 28 

1 0.0 18.878 



APPENDIX B 



SAMPLE INPUT DATA 



o c 



00 

(NivOf^ ^ 
sO^sOO 

(N^O^QO 

C^O^OOCC 

• • t • 



o^O'Ln mC^rooof^sOLnf^'OiAvf'^’n^ovi'— < ir\ \r\ 

m o^G^oo 0^00 c^a^o<7'G^G">^‘C^o^c^oG'C^o^a'a^r^fn 'Oro msor^m 

<rOLno‘^c^ONj-oooo^(NiOLOOLno>rO>j-f^ooa^^rr>aNL'^^r-rnaNLn^fOQOmf^^N^'^N^Ooo 

••in 

r^vC<froi-HO •(^QOscmor^oof^vO^m^o tr^comrvja^r^cD'>o*^Nd'<\i^a'cor^oo«>r>rooooco^-'ONj- 
»-H 1—^ t—t r>-H •— I r-^ CC 00 (-^ •— < «-H i—t f— ^ «—( f— 4 i— I r— 4 «— I 



OOOOOOO COOOOOOOOOOOOO r«-HNtCC^ 

ooooooo ooooooooo^(j'Cr'OooQOoccoQOooaocccoLnLnmLrMrMn>^roro 



OrorncnmmsO'O'OsO'O'OsCOsOooccQOoooo^^— ii-H^^i-HrHLnLrvirMAiAcooocxjQO 

i-H I— I •— ) ^ f-^ *— H t— I «—> «— 4«-H r— 4 I— 1 1 — ( r— I I— » 



<N in 'O 00 O' o ^ (N m Nj- in sO 00 o ^ (\i CO sor^ CO O' o ^ (N ro Nj“ LH sO 00 o ^(N rn > 3 * lA vO 

^ ^ ^ ^ ^ ^ ^ ^ r-4 r-4 CM (\i (M (\i Os, (Nj rg og f\i c\i m m rn fo cA CO ro rri m CO 'T 'J' <“ >4“ Nj- 
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<fo<) 0 rnro^in^r- 4 <j'vcr- 0 v 0 c\joor\jcn 
rOv^> 4 'ininoo(\j^r^orvi>l‘vOOrONi“>rcnLn 
ro ro O' cn f\i nJ- f\i rn c\i vi" ,--1 00 ^ ro c\i r-^ 00 ^ • 
vC >j" ou CO cn >i" 'O 00 CM <M ^ (T^ 00 'O >3" ^ 

oooocor^r^mroo^tno^c^O^oooocoooaoooO^ 
o^c 7 'C>c 7 'C 7 ' 0 ^o'oooor^cr( 0 'c^G'(T'a^o^c>C' 



vJ-cmOo— M"»^ cnmcMM">3'CMOO^ cmm- 
.— <sO^^^^C7'ro^“^Ol/^^C7'c\lO^C^^,— jvOinr--Q'mr^'Om 
• •o^oo(J^oc^^-r*-r-^oao^-ooc^om^cc^Ov^•^^^«r^ • 
^ ^ 00 vO 'O in M- m CMCM CM M- CO cMoo 1 ^ vO nO tn M" m (M r-» 

GOCOoooocoaocooooocCOOOGoaococoaoooaocooo 
C^C7^0^a'a'C^G'C7'O^C7'C^C7'C^C^C^C^C^C^OC^C^C^a^ 



ooo' mtn ^-^ocoN^•— 'r^cMcoroM’LPiC'CMoocM sOooNOr-»vOC 7 'CMinooinrnvCCMoocM*— ‘vOr^r^^inir.cnr^O'O^ 
cr)tncoo^^cocr^ror-roc^f^cM(^ 0 ' 4 ”ror^j oomto*— J oocnro^CT'CMtnf^M’CMcnmrO'Ocr'^cM'OOir^oro 
vO'O'OM-r-HOor^r^oooO'Ooor-ooooo^oooo cMmooi^cnoor-Kfs^sOiroo vOCMocoooor^vOOcnsOOo.— too 



ro CM ^ O 0 ^ 00 nC cn (T' vO <) rn c\i f-t ^ o cn vf CM c nO cn vf ro CM O tn M" CM o r*- nO cn cn M" rn rn CM ^ 

f—< r~t p— ^ f— • 



vOcmoom- cm> 3 -i^co rnc^cn--^r--CMoo> 3 - cm<n^»-^ ^^go cm^^^^OCO u^iTkOtn 'O<^^*-^cmc 0 'O 
cnO^cMvO cMcMCMCMcomo^vCcocT^'OCMCr vOC'O— rnrntnocnvCCM'Oo^ vj" cm^^o ro ro cn m 

(xip-»p-^Q cncnintrur irunvCr^i^ooo^CT' OvOoocT' o— toocT' CT' CT'OOO^CMO Ovo 'O'O MrncMcM— tooh- 



00 00 00 00 00 ( 7 ^ O'O' O' (T'OO oooo OC ^^f-^^CM^<^iCMCMO^mrOrnrOrOM" vf vQvO sO'O'OOO'O OO^ 

p_i ,-1 ^ ^ ^ ^ ^ ^ CM CM CM CM CM CM CM <Ni CM CM CM CM CM CM CM CM CM CM CM CM CM CM CM CM CM CM CM CM CM CM CM CM CM CM CM CM CM CM 



00 0^ o CM M“ in vO 00 CT' o '-t CM CO <r in sO 00 CT' o ^ CM rn vf iTi vO 00 O' o CM cn sf cn vO r*- 00 O' o »--♦ CM CO vt 
vfM“sj*irMninininininininin'ONO'OvOvOvOvC'OONOC*^c^c^i^(^r^r^c^i^i^ooooaooocoooooooooocc 7 '( 7 'C 7 'C'C 7 ' 
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15.949 



^OQOO(^JOC^J<f'0^^'^OCCOCOCrsI^d-c^J^^ncOvOCOOc^JOf^^^^O 

.-H^r\jr\jr0f0rrir0N^>^N^mir\'0ON0i^r^i^f^c0c0(7<7'OOOO 



»^f-^f\i(Nmmmn'i<r tnin'OsO'Or^r^r^i^oocooooooo 



4_H^.-^^^^.-H^^(\jrOs^tnvor^cooo^c\ifON^'iAvor^oo(;ro*-Hr\jfrivrLnvOt^Qo 

^^.-H^(\ic\if\i(NiroN}-Nj-s4-Lmr»mir»NO'0'Or^cocoxcoo^o^O^G^ 



m O' o ^ m (Tunr^ m ^ o •— I m iT\ m iTi r^ CM^ 0 ^ f— j rn 

mo 'Oco ir» f-n c\i c\j m ro m m %t ^ uvr> sC sO x 

m iTi ooi^ 

>j-o^QOo<^.'mimvOaoo»— « 



^ (M ^ sO ir> <r m CM ^ ^ c\i xf vO 00 sC 00 o <\i o (M nJ" o ^ 00 o CO c CM vj- r^j <r vC CO vO 00 o c\i 

t •(MCMrnmmrnsj-Nt >^ir>ir>vO vOvOi^r^r^i^cooooc^ 



,-H ro in uo 0^ r-H 0^ »-H rn un r*'^ in 0^ r^ (7^ < ro ro in un 1-4 

^,-4— 4fM(Mromro>^>^>r'^inm'ONCi^f^f^^ococcoo^ 

m in 

vO'Oinm 



o^oo^ornrnromrororomroo^fMcn^i^Nor^QOO^o^cMrn NMnvOf^Qoo^Of-HnjmNi-msOr-- 
(M CM CM (M m ro m m ro ro ro ro m ^ ^ CM (M CM (M m m nt N^- m in in m nO 'O c) >o CO 00 00 00 cr O' o 
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vOr^oooo 

CT'O'OCT'O 



»-<cMro>^invOr^tnr^o^'-^o— ‘min 
O O O O O O O CM CM CT| m CO 



in CM^ 0^ ^ c^ CO in in 
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r-o 
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^c^in^ 

ooo 



^cMro<rinvDr^oocro»— '(^jrovrinNor^xo^o^cMro<rmsOi^x 
^ ^ ^ r -4 ^ ^ f -4 ^ ,-4 ^ CM CM CM c^j oj n: CM J CM 



rocT'O^'4'm <X) in ro 0^ CO ^ Nj- c\i nO r--i ^ n 1" o CO -4" sO 

Lnh~f\ic>coNOrnr^^co— ‘f^of^a^oof^vocNOoNvco^OsO^ 

lacoOnCi— 

inLni^r-^co^^tLnNOrnmoosj-t^JCNCNJfMogccLn^^cNCNfMr^ 
inLnLn(Njcornr^(^f^c\j0^r-c\jo^r--Lno^rnNj-vor-r-C'Lriroio 
r-- h- vO nO in n 3" r<^ r- vC o nO nO nO (^ h- CO 00 vC 4D 



^^r-^^.--^(^J(^J(^J04^-^♦-^«-Hr-^<r^rnrn^n^•-^.-H*-<rn<^Ii^^cf^ooc^O^^^J^^sj'Lns0^^aoO^O^c^l(^^4•Ln'0^-coo^ 

Ni-^'j-Nj-stvj’s^-tnirunLTLnLnLr.Lr.Lnu'^'OvO^o^^NONOvONONO 



coQ>o^cNrnsj-msOf^coc^o^o4ro>^*Ln4Dr^co 

r— ( »— < 1—4 r— ^ •— < »— < f— t C\J C\j CNJ (^ J C\j C\i (\) C\| f\j 
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'Or^aDCJ'Oro^d-Lr^coC^o--‘(^J^o^d“r^coc^(^J^ONi-ln^o^-cO’-^<^J^^.'0^^ooc^r5^(^Jcn^or^ 

v^">j'vf^^•u^>Lf^LOLf^Lr^Ln^C'sC'Ov0^c^o^o^“^^^^^“^^^“f^ooaocoaococoacc^^<^c^o^o^ 
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NASA TASK-1 TRANSONIC COMPRESSOR 
OF NODES = 107 



APPENDIX C 



SAMPLE OUTPUT LISTING 

(S 

<r 

O 



OCOOOOOOOOOOOOOOCOOOOOOOOOOOOCOOOOOOOOOOOO 

UJOOOOOO OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 
QC 



o 

< 






OOOOOOCOOOOOOC OOOOOOOOOOOOOOOOOOOOOOOOOOOO 

CQOO oo oocooooooooooooooooocooooo oooooooooooo 

< 



00 «-H ♦— I ^ ^ 1—^ t—i *—1 ^ «— ( 1—^ g— ^ f— I f— I f— 1 1— < ^ *— < *— I «— I »— 4 f— 4 *--1 «--4 »— I ♦— 4 *-H *— 4 ^ f— 4 ^ 

01 oo OOOOOOOOOOOOOOOOOOOOOOOOOOOO oooooooooooo 

H“ 

< QOQCiQOOQQOqqCOOOOGGOQOQQQCiQOCQCiCiCiOCiQOQQQQO 
z: oo oooooooooooo OOOOOOOOOOOOOOOOOOOOOOOOOOOO 

—i-^OO oooooooooooo OOOOOOOOOOOOOOOOOOOOOOOOOOOO 
Q.-.00 oooooooooooo OOOOOOOOOOOOOOOOOOOOOOOOOOOO 
Qdwoo OOOOOOOOOOOOOOOOOOOOOOOOOOOO oooooooooooo 

OcQOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 

O »— 4 I »— 4 f»4 f— 4 f— 4 f— 4 ♦— 4 »— 4 f- ^ *--4 f— 4 *— 4 f— 4 f— ^ 4 «--4 fi->l *—4 »— 4 »—4 1— 4 f— 4 ^ f— 4 1— 4 ««>4 »— 4 <—4 *—4 »— 4 ^ *— 4 4 ^ f— 4 1— 4 *—4 f— ^ r— 4 

O • • * • « « • 

OOOOOOOOOOCOOOOOOOOOO OOOOOOOOOOOOOOOOOOOOO 



< 

o 

z: ooooooooooooooooo oooo oooooooooooo OOOOOOOOO 



UL 00000000000000 ) 0000000000000000 . 00000000000)0 
O OOOOOOOOOOOOOOOOOOOOOOOOOOOOCOOOOOCOOOOOOO 

-^coO'CooooLPiO vf oooo cMPiO oooo in oroaDr-ooi^vOLr rooooomoino 

>»-.(^roooooooocoOo 00^000 ooooooNooo'Oina'O^a'O'a'o^nj^or^roOsCo 

Q^w.x>vj-omomoNj-o^>r<^’Oc:'Ovi-rvjo mOir)O<j-a'Nj-r^ 0 Q 0 ^ vrmo oomcomr-m 

<a;cof^NOvfrn '-^0 mo aO'O criOOcof^^O ^fro-^o mocomrvjo ^oDsOtn vj-cvj^O'Ln*-4coin. mOr\j 
CO »*4 ♦—I ^ 1—^ CO »~4 CO *—4 ^ 0^ co co 

OsjX 

O OOOOOOOOOOOOOOOOO 000 o OOOOOOOOO 000 OOOOOOOOO 



II 00 



00 



LU 



UJ 

UJ 



LL 

c 



00000000 00 00 0000 000 000 000000 00 000 



QQQOOOQOQQOQOOOOqOQOOOQQQOQOOQQQC) 
000 OOOOOOOOO 0000 000 00 00 0000 000 000 
o 00 o o o o o o 00 o o o o o o o o <r <r nT vf vT vf 00 o o o 



OCOOOOOO 000 OOOr\Jr\jcvjr\jrNjOO 000 000 000000 
w 000 ooooo 00 00 0000 000 coco CO CO 00 coco coco inmipi mm 

M oooooooooooooooooo^a'^^^.-»^^.- 4 »- 4 ^mm.mmm 

OOOOOOOOO rompncnrONC vO'O'OsOsO^ voo 00 CO CO coco .—I r-j«- 4 ^ ,-4 ,-4 ^,—4 



o c 



00 00 o 00 00 000 00^ 00 000 000000 00 000 00 0000 00 ooc;> 



UJ 

Q^oomvfmsor^cooor- 4 f\jrOvtm>oi^coa' 0 *-^(Mrovtmvor^cooo^r\jm<rinvOf^cooo^(M 
Q ^,--4^,--4^^»--4^^^f\joo(Njr\j(\,(\jnjr\jr\j(NirorOrnrr)rOrnrornrnrnNj-Nf'4" 
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^^OOOOOOO*— '0000.-<00000000 
OOOOOOOOOOOOOOOOOOOOOOO 

oqoqqoooqooooqoooooqooq 

o o ir\ un •-< m CM o m un h- O n m r- 

OOONCor^.— iinroOinoor^c^f^^inOOO'-^cMoorn 

C vO r-* <5 m 0^ CM in 00 o ^ >d- vO > 4 “ CM •— < O CM un 

00 m ro o •-< vO ro o r- •-< r- o 00 un CT' ^ m 

ooo^u^.oooiT\cMO^Ovomor^or->tvMnLrMnso<5 

ia'C'aNcocoQcr^*-<o^oor-'vf*-<(7'a^oof^vOLnm»-<ooooocooooooooooooooooooo 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 



OOOOoOOOOOOOOOOOOOO OOOOOOOOOOOOOOOOOOOOOOO 
OOOOOOOOOOOOOOOOOOOOOOOOOOOOOCOOOOOOOOOOOO 

CtOOOOOOOOOQQQQCQOQCiCiQOOOQCiCiQQOOOQOQCiCQQOOO 
f^oON^rnof^h-f^oor^t^coc^co— <>d“(MvO*-^»-^<fvOoo>tvOOr-oo^Nd-cMvCvr'^'^f^ocMoor^ 
inr-cMOoovOc<^h-^oo^t^Of^croor-vocMOOvooOsO’“‘vOOOOf^vC<fcMM“vo<7'Or^ooo 
LA m o o vM^ 00 ^ o >t ^ m <f ro o tn o o ^ nJ- ^ CM CT' r<^ nJ- ro tn o >t (T' o CM o o 

Lnunm.— 100 '-^'J'LH 'OmmOO'4"CMCMCMCMCMCOm^*-^CMCMCM(^'-<CMCMCMCMCMOOOO.— isCvOvOC^in 

LH CM 00 rn rvj o •--< CM (T' r- Ln o ro no o LH rn LPi r-- o LH o >t CM >o o ro 

onr-r^i^vO'OvO vO'O vof^r^r^ oocc vO'Or^r^ oono vOvOsO^^r-r^ oooo rnroco^tsj-cooooo 

OOOOOOOOOCOOOOOOOOOOOOOOOOOOOOOOOOOOOOO ooooocooo 



OOOOO ooo oooo o oooo oooo oo—'^^^^o OOOOOOOOOOOOOOOOOOO 
O oooo oooo oooo oooo oooo oooo oooo OOOOOOOOOOOOOOOOOOO 

QOOQQQOQQQQOOqQQQQQOOOQOQOOOOqOOQOQOQQQQQQQCiQQOC 
ro sD ^ o vO o m m ^ m ^ vT/f^ o vo cMoo cMin oo o c o vi" cMO o ^ vj-r- in si- mcM >d- Nt CM o o ^ o 

vO *-<vOOm vMnLnoocM*-Hr- Ocm nJ- voom >d“vf inmoo ooo^ NOm,r^omr-sOin'-<o cMOo^vOmf^ocn 
•—<0 00 *-<00 mmOLncM>d“ cMcn cM>^* —^00 *-<mtM^ooir. oooooomvooocT^r-r^r^vOoof^cooomvOOOM" 
cM^OoovO n^cmoo invf •-<nCoocmcm ^ooO'ONt(MOr-.ooooocor-r^vOvOinvfro(MCMCM'4"<^.cMcoh-h-vOvOin 
o ooo oc-oo oocor^f^inm oinooNCJ^ oooocoooco coo oo ooooooococooococooocqooooooococooooooo 
OO 0^OO(7'OC^OC7'OC0C0C^O^O0^OO(3^OO— ^oc^c^c^ooc^c^c^c^oc^c^c^c^c^oc^c^c^‘ 



OOOOOOOOO oooo oooo oooo oooo oooo OOOOOOOOOOOOOOOOOOO 



(MCMCMCM(MCMCMCM^-<CMCMr\JCM'-<CMCMCM(M(M(MCMCM(M CMCMCMCMCMCMCMCMCMCMCMCMCMCMCMcMCMCMCMCMCMCMCMCMCM 
OOOOO OOO OOOO OOOO OOOOOOOOO OOOOOOOOOOOOOOOOOOO OOOO 

Q OOQQ QOO O oooo OOQOCQOOQOCQQOO oooo oooo OQQCiQQCO coo o 
OOO OOOOO OOOOO OOOO OOOO OOOO OOO OOOOO OOO OOO OOOOOOOOO 
o^ CM in 00 oomoOf^vCooo^h-cMoocn ^tmOcMOOcMO'^co vO'-^vOOcMLncomro vOcmoo cM^vOr^r- — <imn 
invOf^‘^cnLn(X)Oin<-Hcnof^'4“f^fnocof^cMf^o>j-rncMooo<frn^coincni--iOcMLncn>j“CMinrnrnsCC>*-<cMNO 
(M^Oco^OvOvo<^•cMOo^-'^^mcoC^ooo^^coaDO^coooocMfnoo^-^ncoF-^^tvOvOtnoo^ocMoococo^-.vOOmo 
oor^ s0^rocM.-<o »-<00'Om vOnO^cocm^ ^or^Ln<^cMo^-OLn^4" cm or- m^fcMOr^ sn mm 

r— ) F—l *— I < I— 1 1 — I *— < J F— < r— ) •— I r— < < i— I r— < f— I *— < f— < f-H •— < i— I f— < f— I i— < f— 4 i— < i— I f— < f— < <— * f— < •— 1 1— 1 1— I r-H 



OOOOO OOO oooo OOOOO OOOO oooo OOO OOO OOOOO OOOO OOOO oooo 



(M(MCM(M(MCM(MCMfMCMfMCMCMCMCMfM<MC^JCMCMCMCM(MCMCMCM(MCMCMCMCMCM(MCMCM(MCMCMCM(M(MCMCMf\jCMCM(MCM 
OOOOOOOOOOOOOOOOOOOOOOO OOO OOO OOOOOOOOOOOOOOOOOOO 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOQOOOOOOOOOOOOOO 
oooo OOOO OOOO OOO OOO oooo oooo OOO OOO oooo OOO OOO OOO OOO 
^r-ma"OcMco<focMvfr^coomcrin^f^cMoo>tocMcM^— *oo>toocMoocMsOccoinmOLnoOvomcM»-^cM 
^>d“00 F-<maNCMv0 ocMCMcMCMcnmo vomcrvOCMo soo OF-^cnmmom sOcMvOO sj-ocM*-<'-<ooornr-mcrm 
m^4"romcM*— <•— <oommLnmmmmsO^-^-ccc7^(7^ OvOcoo^Of-hcoctoo^ coo^cmvOvOvOO'-On^cocmcm^^ 
<X) o CO oD CO 00 CO 00 00 o o cr cr o o oo o o o o F--< ^ r-H CMCM CM CM CM CM rn rn rn rn rn vT sj- vT o 'O vO nO nO vC 

^ ^ F-< ^ f-h ^ F-< f-h ^ ^ .-I ^ ,--4 —I rj CM CM CM CM CM CM CM CM CM CM CM CM CM CM CM CM CM CM CM CM CM CM CM CM CM CM CM CM CM CM CM CM CM 



OOOOOOOOOOOOOOOOOOO OOO OOOOO OOO COO OOOO OOO OOOOO OOOO 



cnNrmsO!^o<^o,--<cMrnvfmvor^coc>o^cMcnNj-invor-cDOO'— <cMroM-invor-coC7'o*— <cMfnsj-invor^ooc^o 
<r>d“>i'vrNt'4'>tmmmLrtmminmmLnvO>C'C'OvO sC'O'OvO ^or-r^r-r-r-r-r-r-r-r-oo co ex; xcc oooo o coco cT' 



108 



IL 

ocooooooooooooooo a 

OOOOOOOOOOOOOOOOO OJ 

Q- 



OOOOOOOOOOOOOOOOO 



sOQCO(^lOC\|Nt' 0 <^^OQOOOOOC^J<^C^JN 3 -vOoO^OCOO^^JOc^J■ 4 •vO 

^^{N(\l^o^n^^m>^xJ•^J-lJnlJn^o^O' 0 ^-^-^-^-cococ^C^oooo 



OOOOOOOOOOOOOOOOO 

— 4 --4 (M (M ro ro rr, rn ^ vj- mm sO ^ o r- r- r-- CO oc o o o o o o 

»-- 4 «— < 1—4 fmm4 

C*~'j 1—4 *— 4 ,aa 4 f — 4 i — 4 ,_b 4 i — 4 j — 4 f— J i— 4 «— 4 *— 4 

OOOOOOOOOOOOOOOOO 



qoooocqqoqoooqoqo 
ocN<roooooooooooooo 
r-vomoooooooooooooo 
r-r-r-oooooooooooooo 
sT f<^. CMO o o o o o c o o o o o o o 
000000 oooooooooooooo 

QNCT'O^r— I— 4 «— tf— »*—'*-- 4 r-(i-Hr-lf— (f— If— ».— 4 ^ 



OOOOOOOOOOOOOOOOO 



^c\jro>j“mvOi^aDOO«-H(\jm<fm.vOr^cooO’-^cNfnN^'mvor^co 

^ ^ ^ ^ (\i eg cvj rn nT vt m m m m vO ^ nO CO 00 00 CO o c^ o o 



rome-o^h*o»— 'm*-Hromr^mr^c^»— lO^mmcomi^Oh-O^ro 
r- 4 ^(\jcsjrnrornro>^'>^*‘g"mmsO '0 0 ^-h-i^r^oc!aoG^O 



>- 

o 

c\jcNc\jrgcsjr\j{NCNcNj(M{NjcNrv!(\){NjCMc\jO 

OOOOOOOOOOOOOOOOO-J 

o 

OOOOCOOQOoOCiOOOCOCL 

oooooooooooooooooo 
mr^ooo^moOsOooooooomoK 
omoco<f’-^ONj-mrHomOmooor- 
o o ^00 c^ nJ" o ^ 00 o (N m m vc CO 1 21 
coegr-^e-mv^-cNj— ^r^i^s 0 m' 4 ‘mrvj^r-^uj 

f — 4 f — 4 1—4 r— ^ 1—4 *— 4 f — 4 »— 4 i — 4 i — 4 «— 4 r »4 f— ^ 



<Msj-sOCOvOOOO{NOCN>^NOsj*vOCOOCOO(MN3-r\jN3-vOCOsOCOOC\l 
.-^.- 4 c\j(Mrnmrorn>r>^>tmm. o>ovor^r-r-r-*oooooo 



^ m m m e - »-4 m m ro m o ^ m ^ ro m m. 0^ •— I 

r-ir-i»-'(NjcNJrnmrOvf sj- <f> 4 -mm vovoi^r-f^r^cooocoo 



?•••••••••••••«• *00 (>0 

ooooooooooo 000000> UJ 

00 cOr-icNjmvf mgor^cooo^<MroNrm'Or^coo^o»- 4 (\im>i'm 

C^,- 4 ^.- 4 (M(NjcM{\i(^m> 4 -Nj-*nmm.m'ONOvogDoocoQOcooc^o>o 

(N eg (N j <NI C\J CNJ <\l C\l eg (N <N 1 og <\l (XI CNJ CNJ (Nl 
OOOOOOOOOOOOOOOOO 



QOCOOQOOQOOOOOOQO 
OOOOOOOOOOOOOOOOO 
OOOOOOOOOOOOOOOOO 
^moomou^oooooooooo 
r-ioor-vovommoooocoooo 
g 3 sO vOO OO OO corn corn rororornro 
eg (Nj cN (N eg eg CN eg ro ro ro rn ro ro on rn rn 



me-O’— ^Of— ife>mrome-or*»o^ro»-- 4 fome*me-o^f— lo^rom 
^^^egcgmmrn'^'>^''«t>^'mmgDgDe-e-e“e-ooooooa'a '000 



OOOOOOOOOOOOOOOOO 



Lu *-4 eg ro m ' o 00 o o eg ro m sO h- 00 0^ o — 4 eg CO m sO CO 
21 ^.- 4 ^,- 4 .- 4 f-^r- 4 r-H^^egegegegegegegegeg 



eg ro vt m vO 00 cr o eg ro vt m go e^ 

c^a^oc^c^c^c^c^c^oooooooo 



UJ 



00 





























o 
















































CM 


























0 




o 




















0 






LL 




















0 


a 






























oo 




















Q 


tn 




< 




















0 


CM 
























0 


r ~4 




LU 




















0 






a: 








LU 












0 


r -4 




< 








X 












>t 


• 
























nj 


0 




uo 














0 






f r— 4 






LU 








X 




0 


LU 




r -4 


00 






_J 








< 




LU 


C/> 




0 


1 


II 




CD 




K 




q: 




LO 


\ 






Q 






< 




LL 






X 


\ 


I— 




Q 


11 0 


X 












0 


CL 


h- 


LL 




0 


0 


a 




a: 




3 




LU 


QC 


LL 






0 


LUO 


*— « 




< 


0 


0 




Q 






ro 




0 


oca 


f— 




> 


LU 




H -4 




<t 


ro 


0 




0 


xo 


X 








X 


on 


ro 


0 


C 








00 00 


X 








CD 


a 


0 






Q 




4 


00 * 


H -4 






X 


_1 






Q 


C 


1 ^ 




• 


LUO 


QLQ 






CD 




CN 


0 


0 


CO 


vd- 


(M 


0 


cc 


LU 


(MCMCMCMCMCMOJCMCMCMCMCMCM 


< 


-4 


f-H 


0 


0 


m 


OD 


CD 


0 




Q- 


00 »— 1 


ooooooooooooo 


2: 




0 




OD 


0^ 




CO 






11 


LL 




> 


ro 


1 


Q 


ro 


m 


1 -^ 


1—4 


Q 


11 


1— 




QQOCQQQQOQQCC 


Q 


0 


Q 




O' 


ro 


0 


ro 


0 




Xf- 


0 


«— 4 1-H f —4 r — 1 ^ *«-» 1 -^ 1-H 4—4 


0 




0 


VsO 


O' 


<r 


rvj 


• 


0 


00 


<r u i 


►-^LU 


♦-MLr.tx^tnLnLntnLncninLnin-uoLo 


2: 


0 


0 


< 3 - 




• 


• 


C 


LO 


H 




00 0 - 


^ CM CM CM (M CMCMCMCMCMCMCMCMCM 


q: 


0 


CM 


LO 


• 


0 


0 




ro 


< 




aoo 


4 4 *^ «— H •— H r**H 1— ^ 


LU 


0 


0 


0 

in 

1— ^ 


0 








ro 


LU 


X .-4 






X 


0 










II 


in 


X 


c 


LLOO 


Q_ 4 ^ f —4 »-H f —4 4 -^ 1— H 4— ^ f —4 1— H 4—4 f —4 4—4 O O O O O O O 




f^ 


CD 


• 




II 


II 




• 




o^- 








0 


• 


c 


II 






>- 


0 


0 


<r 




ooooooooooooooooooco 


1— 




C 






Q 


> 


1— 




1—4 


h- 


LUt -4 




LU 


• 






LU 


LU 


f— 


*— t 




LL 


<i> 


f-oo 




—i 


0 




II 


cc 


LU 


♦-M 


0 


II 


>—4 


H- 


<r 0- 








II 




X 


Q- 


0 


0 




0 


1—^ 


X 










LU 


K 


LT) 


0 


_J 


h- 


LU 


<00 


»-hLU 






ii 


> 


a: 


< 




_J 


LU 


X 


Q- 


LU 2 : 


K q: 








h- 


X 


cc 


-J 


LU 


> 


< 


00 


XLU 


00 LU 






LU 




00 


LU 


< 


> 




F— 




Q 


LUX 








CO 


UO 


CL 


X 




X 


00 


LL 


0 


3: 


• 




< 


21 


LU 


X 


a 


X 




X 


O 




—1 


0 




a: 


LU 


q: 


LU 


*— 1 




I— 


0 




LL*-i 


<00 


X 






Q 


a 


1— 


H- 


K 


LU 


0 


0 


*- 4 |— 


»— 4 LU 


0 in nT O' 00 ro (M 1 ^ 0 4 -^ 0 in <r ro CO r- CM 1 vO 




3 








< 


LU 


_J 




»— 4 


o< 


h-O 


yj,— (CMcMroM-LnLT or*-ooooC 7 ''^cMCMro> 4 'LnLn 




Q 


t— 


f— 


h- 


1 - 


_l 


I— 




1— 


LU f— 


-mO 


0 




—1 


0 


0 


0 


0 


X 


X 


< 


< 


Q-i/) 


XX 


g 




LL 


K- 


h- 


H- 


QC 


»— * 


c 


0 


CC 


00 


»— 4 


X 



110 



I— < o 

(NirONj-or^oDO^ 

o 

o 

Q 

o 

o 

o fo sO o (Ni rn vO r^ 

(Nm«4-inr^coa' CO 

r^ 

rn 

m 



c^ 

c\j 

o (\i m CO 0^ (\i LH >j- 

^rOvTinsOOOO >j- 

o 



c\j O 

o On 



J (*Vj O' J T— ^ T— ^ T— ^ f— ^ OJ OJ OJ 

oooooooo oooooooo 



QCOQOQQO 
^ OJ ^ OJ o- (\i >o 

u^o^-s0<>’-'(\icn 
ojooor^^voo 
^ojinvrrooor^r^ 

r^^^o^vOtnr^csj 

OOoOOO^’-^t-^CO'OnJ-oj^C 



CCQOQOCO 

^0^rO'4“rO— 

LTiOm^or^>j*Nj- 

OjOOC 0O--^O^ 

^ r-H *-< o- 'r\ o o- 

r^vj-ojcoLnNj-NTvO 



LL 

o 

tu 

a. 

oo 

U1 



oooooooooooooooocooooooo •- 



LL 

LU 

CL 

LU 

X 

3 



in o vt rn GO r-i OJ m LPi nO CO (7^ o o '-H CM rn LH vc LU 

sor^r^cooo a^ooooooooo 

(— » 1—4 1—4 r-^ <— J ^—4 •— 4 O 

z: 



00 t- 4 ^ un CO (\i O 


Q X 


^m>j-LnvOCcc^o 


0 O 


OJ 


* 0 «J 


in 


0 KH 




LU 00 00 


o 


•— O. Q_ 


CO 


LL rO LU 




00 rn 




^ LU 


^rO'^in'Of^OO 


)— 0 0 


sO 


< csjoj;z 


m 


OO nJ- IL* 




0 


• 


• • qC 


o 


LU OC'JJ 




> :z> 


vOo-Om>j-r^o 


X 


^oj-<finv0r^a^oo 


K ooXO 


»— ■ 


c <— Oo 




^00 


r— ^ 


OOojf— CL 




lulu < 0 




LuOa: q:ll 


2I-Z02UJ 


rOOO'OrOvOO'CUiLUCOt— UJ 


^ojminNor^co»— 




HQ-CLQlH- <l 


<3:01 <z: 


cLO>:zcLC:2 


luo^>-*lu n 


K 


0 h~Q-i 


t—* 


I LU _j 


OJ LPi o C> OJ in QC 


hh i-c 


T—tojrnvfOo-oocL 


^ CL<LL 


o 


ozrox 


LL 


• ^CLL-^X 




0 1 — 2 : L_; 


00 :z OH- 1 /^ 0 ^*— 


a. 


^<ClluX 


Lu:23crLLK-3 


^ OJ in CO ^ 


OlLLU 


^OJCOnJ-sUO-COH- 


< h— 1— ‘ 2. 00 


OO 2 : >-< 00 < h“ 


LL. < < LU QL -J 


OCdUJK-OOX 


C^LUQCXCCOOO 


<J— H LU< QLLL 


— J 


>-^oo^^CLa: 



111 



FINITE ELEMENT RESULTS 



OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 
I I I I I I I I I ) I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I r I I 1 I 
>-QOOCOOQQOOOOCiQCQCiCiOQOOCiOOC;OCiOQCOQOQCiOCOOCQQQQ 

I— o o o o o o oo o r^oD m nO inm vo nO nO oD rvjr^ m vf CO ct'cnj vD O' ^o o rvi vo 'O 00 ^ cMCT LA o 

o oo ooo ooorna'cor>-aoOf\icMO^Or-ooa' sOooroinsOiA'O'— «or^ 'O'OmcMp-Hor^G^’^^cMu^ 
L>oo oo oooocor-(\j vCocoo^ ^^crmo^pvOCMiA'OiAcr vOr->tro>t'OO^fM< 7 ^>toa'omr^f^o^ 
-iiooo ooo oooc\iv 4 “>tinNtmfornroro^LA>r vf •— Aj^in^miA'-^oo (Mr^cMr^omciAC'OOin 
Luooooooooooooooooooooooooooooooa'Oooa'crOo^oaor-{^rn>i"m 
OcoooaoaoaoaoaoaDooaoaoaocoaOaoaoaocoaococooCcoaoaoaoaoaoaocof^r*~r^r^r^aoaocoaof^r^r^oocooo 



oooooooooooooooooooopooooocoooooooooooooooooo 



<^J C\J (^J (\J C\J C\J (\| ^ ^ C\J C\J AJ c\j ^ (\j <^J (\J <\J (^J CNJ (\J ^ C\J f\l (Nj ^ (NJ AJ (\| (^J OvJ CNJ ^ 1-^ ^ C\l (NJ C\J CNJ ^ (\l C\j (\| 
OOOOOOOOOOOOOOOOOOOOOOOOOOOCOOOOOOOOOOOOOCOOO 

QQOQQQQQQQQOQQQCQOQOCiQCQQOCiQQQQQQCCiCiCOQQOOClQQ 
OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO o 

-^co (T' o o oo oiAO o c oo O' mo o oo oiA o cA 00 A-o or- o m vf m o o oo LA CLA o o 
r- fA o o o o o cr cr CO o o o cr o o o o o o o CT' O' o O' cr >o LA P cr p P O' Aj ^ o ro o o o LA o A- 
'^oOvrOLAOLAOvTcr>j" 00 ocT'NtAjOLAOLAo^o^vfr-ao(T'vrmpLA^r-rACro^a'maoror-rAAj^o 
Qcaor-sO<frA»— ‘OLAOcoNomoc cor-sO vffA^OLAoooLAAjo^^co nOla vrf\ii-HO'iAf-<aoLArnoc\icor-NO 
^ *-H t-^ 1—i r— ( 1-^ f— ^ 00 r— •— I f— < < r— 1-^ »— 1 1— < p— ^ 1-^ 1-^ IP — 4 CO f*** i-h f — i i— h r— 1 1— ^ i— < i-^ i-^ r— ^ (y^ oo r~ i— ^ i — h ^ oo r-^ 



OOOOOOOOOOOOOOOOOCOCOOOOOOOOOOOOOOOOOOOCOOOOO 



ac 

> 



ooooooooo 



ocoo OOOOOOOOOOOOOOOOOOOO 



oqqo 

rom LAO 
^r-LArn 
O O *— • rA 

AJvTO AJ 

r- «— • cA r^ 



oqoo 

vDLA^rA 
^^00 CO 
1 — (•— if^OO 
COLT* 00 O' 
CAsOLAA- 



^r-AjOvf>i’»— 'vO 



QQ 

vOsO 

AJO 

>J“rA 

cor- 

sOsO 



QCQOQQOQQQCiOOO 

«r-cococroAj'-<Aj»-HiAiAONO 

A J A* LA AJ CM la rA AJ CA LA LA AJ nO AJ 

r-cMuAr-o-— •vDOfA'Ox^r^vCoo 
vfLAsOONOaDrAocrLAsOAjcro 
vOAJrAr-sOvTAj^ocrsOA-r-rA 
LACMrA'J“^^AJfAi— ir-^^ 1 — 



CM CM Asl *— i CM AJ CA CM CM CM 

COOOOOOOCOO 

I 

OQOOQOQOqOQ 
COlAOOO^OOOCOCT'A-vDa' 
O »-H LA AJ O A- CM AJCO ^ 

lArAr^OiANj-O^Ocop-o 
cMrAO^O'CO>d“a'>fCMrACO 
nO CO sC O CO AJ ^ ^ vT >0 LT 
vO^>i"AJvf^>t^LArACA 



ooooooooooo ooooooooo ooooooooooooo ooo ooooooooo 

I i I I I I I I I I lilt I I I I I I I I ill 



rArArArArArArACArAirArArArArArArArArArArAfArArArArArArArArArAfArArArArArArACMrACArArArArArA 

coo OOO OOO ooo ooooooooooo ooooooooooo ooooooooooo 



QOQQQQ 
cArArArArArA 
00 OD 00 CO 00 00 
fvj>t Ni"vT vT si“vT 

y> 1 -^ I— 1 1 -^ I— I f— ^ I— j 

oooooo 

CM CM CM CM CM CM 



COOOCQQC)CCjCiQQQOCjOQCCQQcAC.OQQOOQQCQCjCQQOQ 
CA rA rA 0"0 CM •— « A- CM CA 'A\ LA A- rA o ^ la CO rA sO o cr O O r* 00 CA vO CA O CM A- LA CA r- 

coaoocLArncMCMsj-or--’ 4 -aO'Oco>tO>d“r-aOsOCrOAj.-irAoocoLAr-r-^^coCMLAcM»-<ooA- 
vr^r'tOO^'OrAPLALALACOCMvOCrO'OsOCMrArAvfsOLA^OOsOaOCMfAA-CMrA^^a'O'rA 
^^^OONDCMLA^sONOAJf<">CMNfvDOA-LAO^OCO' 0 <^^CO^LALAOO^rA^^CrLAOOCM^ 
O O O A i — < O O O ^ ^ ^ O O ^ O CM AJ (JV AJ GO O AJ vf vf N^* I— 1 00 LA LA sO 00 CO OD 

CMCMCMCMCMCMCMCMCMCMCMCMAJCMCMCMCMfMCMCMCM^CM^AJCMCMCMCMCM^O'AjrArArA^f^f-^ 



OOOOOOOOO OOOOOOOOOOOOO OOO OOOOOOOOOOOOOOOOOOOO 



CM CM CM i ^ CMCM«— 

OOOOOoOO oooo 



CM CM CM I— 1 1 — < I— 4 1—4 

OOOOOOOc 



CMCM^^ CMCMCM^^^^^ CMCMAJ^ AJCMCM 

OOOO OOOOOOOO OOOO OOO 



^ooooqqqq oo 

^ I-H CM ^ CM r- r- CM nO •—‘LA 

-^lAOA- vCO^CMCA LA^ 
^Ajoo 'Or-p-^sO c Aj^t 
cO_iCMLA>J-rA00r-A« ^CO 
Q_(«^ vf ^crvOinr-cM r-*-^ 
4 »— ii— lOO'OvfCMi— 



QQ QOQQQOQQ 
sOrA ^OOO A-<MCMA-0 
LAO LACAOOOD AJIALAO 
CMLA AJvO^CO voor-o 

OOAJ •^LA,-<r^LAr-v 0 CM 

A-00 r- CM cA p r- 00 rn 

sOAJO^^^CTvD>J“AJ,— I 



O 



QOOQ OOQOQQQQ 

•^LArAvT 1 — I sO^LAlA't vTOO 
LArACAvO LAA->f O-^CM A- vT 



cMr— *^'0 CM r— cA ^ 00 vT 

^O'CMAJ CM CM <r CO CM vf 

r-^oo' r-vTAjr-Ajcro^cM 

^ , r- Aj o '-I ^ cr A- vj- Aj ^ o 



oooo QQQ 

^XA-^ ^XCM 
LA.-«>tO LACACM 
CMf-^NfvO CMCMCM 

A-LAOvO A-LArA 

^ Vj* 1 “^ 



coo ooo ooo ooo 0 0^00 <^000 0000 OOOO OOOO oooo OOOO oooo 



iU 

AJrA >tLAvOr^CCC 7 " O^CMrA vTlAvO A-COCrO^CM-'A vTLA'OA-XO' Oi-<AjrAvrLANOr-Xa^Or- 4 AirA>riA 
O ^^^^^^^p.^r- 4 ^CMCMA^lCMCMCMCMAJCMCMfA^rArA^AcA^A^A^A^A^AN^s^•v^v^v^Nf 

2 ' 



112 



OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 
I I I I I I I I I I I I I I I I I ^ I I I I I I I I I I i I I I I I t I I I i M I I I I I I I 

CQQQQQQCOOQQQOQOOCQQQQQQCQOQCiOQQCQCOQQCQC-OQQQcQQ 
Nj-<Munr^oo (Nj(\j^^ooo«-^ma'a>NOO ^rvjmcocoo^ vO 

oDaDi^cn ' 4 “U^Or-^-o^a' 0 ^-a^irun'Omino^aDOOOo< 7 '^-^-o^sOrMrnaD^<\iaDO*' vf^r^rvj^-uMno^NOONrr^ 
<) sT ^ f\i o ^ rn ro c in rn sj- ro 00 ro ro tr h- vO tn o O O ^ tn nO o ^ ro (\j ro o oD ^ ur O ^ 

vj- ^ ^00 o ^ o vO vO o o o m ro r-^ 00 o m (Mm (Mh- r-H n 4 “ m (MO 00 o mf-H o^ o m vC o m m m ^ m vO 

o r- m ^ 00 o (M (M o o o — < CM 00 ^ 00 o <5 oD vO Nj- vO vO r- o i^^' m m (M o h- vO m cc oD m vt m m 

oD so 00 00 1 ^ vO <) 00 CO 00 vO vO vO nO h- h- kO 



OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 



(MCM(M(M(M^CM(M(NJ(M^(MfMnJ(M(MCM(MCMfM(Mnj(Mnj(M(M(MCMfM(MCM(M(M(MOsJ(M(MOsJ(M(MCM(M(M(MCM(MCM(M 

ooooooooocooooooooo OOOOOOOOCOOOOOOOOOOOOOOOOOOOO 

qqqqccqqqcqqqocooqcqqqqqoqqqoqqcoqqcqqqqqqqqccqc 

OOOOOOOOOOOOCOOOCOCOOOOCCOOOOOOOCOOOOOOOOOOOOCOC 
m CD G^ o ro o o vO 00 O ^ h- (M CO sj- m O' (M 00 nj o <5 00 vO 0^ nj m 00 m ro vO (M 00 ^ vC r- ^ m m m G^ 

mmmoo om^rnof^'tr^mcoo(^(Mj^o sT ro(Mooo '^fT, ^oDmrOi^G'Ov;mmNj-cMm(^^rovOG'^(MNOOmo 
oovOvOv 0 ^o^oor^r^mooG^oooi^oocoo^ooooo(Mmoor^moo^> 4 "vGvomooNO(Moooooo^-voOroooo^ 
Ni-rncM^o^cDOro^NOr^vOvONj-ro(M^^Op^msj-(MOr-vOm''t> 4 "rn(M^or^mNj-(MOi^sommNi-roro(M^ 
f— ^ G' G^ ***^ 



OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 



fM(MCM^(MrO(\j^CM rorn(MCMfM(MrOrOrnrom^(M(M(M(M(M<M<M(M(M(M(M(M(M (M (M(M CM O rO rorO f^. rO ro rO 

OOCOCOOOCOOOOOCOOOOOOOOOC^O OCOOOOO ooocoooococoo 



CQOCC;aQOC.OQQOQOQQOQ 
^ ro sO m (Mm ro sj- o rno m st P^ G' 

oxG'P^mcMm^i^O ^(MCM«-< vfm 

(MP^r^(Mf^(Mvj-G'NO^^sOOGrnOmmNd" 
(M^OO^O^vj-st(^r^f^P^h-^G'(M- 4 '<Nlor- 
ro (N (M o CM m sj- r- o (M ^ fo (M — I m m 
fO fO «— I CO 'O ^ m <f" (M *— < m oo ^ «— < 



CQQQQQQQOOtHOQQQ 
>d- moo ^ o nO on ^m (MOO 00 
ONmooooooor^r^rOvomOmm 
>d"G'rOrn'-^>J-(<'.ror^r^ 00 (MCM(M(M 
G'OG'(MO'Or^mcDCMG'CM> 4 " 00 P^ 
omvrooomoooommrocMmG' 

^ >r moo G^ 00 ro ^ CM vT m vO P^ p^ p^ o 



qqcqqocqqoooq 

v^sj-r 0 ^ 0 'LPk>^ 0000 vfOC'j 00 

^ (M vC »-< nM^ cM^ o 00 ^ (^ 
OP^r^omc^OOCO< 7 ^CM^^OOO 
r>- sj- > 4 “ vO vO p^ CM nO 00 
vOrOf^P^'GcMP^rnc'm.-^'OCM 
(M>^'in(^m^(Mromrom^ro 



OOOOOOOOOOOCOOOOOOOOOOCOOOOOOOOOOOOOOOOOOOOOOOOO 
III I I I I I I I I I I 



rn rn m rn m m rO rO rn rO rO rO ro ro ro rn rO (^ ro rO (tl rO rO rO m fO m m ro ro rO ro m ro rO (tn ro m ro m rO ro ro ro rO ro rO 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOCOOOOOOOOOOOOOOOOO 



CQC( 3 QCC,QOCC^C:QQQQQQOQOQCQQc=iQOC:QC;COQ(Z:QOCcD 
00 CM 00 00 O m >j- m m ^ ^ rn sj- P^ r- ^ vj- o vj- nO ^ nO m 00 G' <n ro h- sO G' nO 1 ^ CM 00 ro G' P^ 

o o CM m sf o rn Nj- 00 m ^ (mp^ G' ro G' (M nO m o vo '-^m P^ m o r -4 sO G' sO 00 00 p^ o O' 
m o G' vO ^ ^ Nj- 0^ m CM m rn 1^ o 00 00 ro m m vO m < 7 > CO 00 CD 00 (M sj- (M vO o 

f^rnOf^mm(MOCNtoOsj-^NOp^(Mro^,-^-Hp^(MmOvj-Ln^cMG'OoovOoof^.-^oor-cof^o 
nO o m 00 rn vj- O m 00 00 ro ^ m c sj- o (7> ^o o m o ro st- o (G' ( 7 ^ ^ (M (M (M m o m m m G' 
(M (M ro ro vT CM (M m m sO CM CM CM CM ro rn vf vj- IT. m CM (^ m m (M CM ro rn ro m ro vj" m P^ m m 



QQOCOC:iCOC 

G' 00 O(Mr 0 G'mO 00 

OCMOG'P^P^f^rOG' 

'(no'Gmcovj-m 

(MNomsj-(MCG'rnr^ 

^f-^rO(M(M(M^rom 

rO/r^rOrOrOPOrorOrO 



OOOCOOOOOOO OOOOOOOOOOOC OOOOOOCOOOOOOOOOOOOCOOOOO 



(M • 




CMCM^^ 


(MCM(MCM^^^^ 


(M(M^^ 


CMCM(M*^*— 


(MCM^»-^ 


CM CM 


(M 1— H 


oooo 


o 


oooo 


OOOOCOOO 


Oooo 


oooooo oo 


oooo 


oo 


oooooo 


QQQ(=j 


Q 


QOQQ 


QQQOOQCQ 


QQOQ 


qOqqqqqq 


otnac: 


QC 


QQGQQQ 


GO(M^ 


G 


^Nj-P^O 


sj- sj- m ro m 


•— <«— lOP^ 


'-^CMPOOOOvj-P^ 


— 'i^^m 


^CM 


cMm Nf Nj-moo 


CD •— <P^ >4" 


00 


mrOsj-O 


mNOOCM^P^CMOO 


moomcM 


m om ro n4“(M G nO 


mmvO^ 


moo 


OgG(X^O o 


o^mGst 


CM 


CMCMGP^ 


cMf^ G cMm o CMro 


CMnCOO—* 


CMf^DmvONt— ‘O 


(MsO oo 


CMO 


--^roO^roCM 


GM“<fO 


r- 


^CMro O 


'-^P^OOsOf^GOOf^ 


^ooor- 


^OCMOM-'OOOO 


--•OCMO 


•—< 0 ^ 


G (MP^P^sj-rO 


OP^COG 


m 


p^mcop^ 


r-Nf(MOCM^OO 


P^(MG'-' 


p^sj-^mp^Ovj-(M 


P^^>4"CM 


P^ro 


OCCMO'-^O 



^ 00 vC ^ CM O ^ 00 M" O ^ ^ ^ ^ 00 ^ CM O ^ ^ P^ <• O ^ ^ ^ G' P^ vO <f (M O P^ nJ* O ^ ^ ^ G' m sj- (M C 



OOOI-^OOOOOOOOOOOOCOOOOOOOOOOOOOCCOOOOOOOOOOOOOOOO 



vCf^coGo^(MmNrm^oi^coGC->^cM{^.<rir. or^ooGO '-^(Mpo vrm<?i^ooG o^cMn-^ vj-mvOP^ooG o—^njro 
sf ^vj- vj'iTvm mminmmmmm vOO nO vO o vO^CvC'O^OP^^^^-P^P^P^P^p^P^^*aocooooooO(DOOooooccGGGc^ 



113 



f 



OOOOOOOOCOOOOO 
I i I I I I I I I I I I I I 
QOQQQCOOOOOOCiO 

sOOCNJoomrncNJinf^^^otn 
LTicooOsOroh-vO>ff^' 0 >trnm 
•— •OiT\oCsOcx:h-ir\>^’{\jrOs^r*^co 
cor^f^vONCf^r^f^r^f^^*r^r*-p^ 



oooooooooooooo 



<NJ (NJ (NJ CNJ (NJ CNJ CN (NJ (\J (Nj <NJ <\J C\J 
OOOOOOOOCOOOOO 

QOOOCOQQOOCiQOQ 

OOOOOOOOOOOOOO 

oornoosooooooootno 
^ 0 <^—'O'«^’m'-^OlnOmocD^- 
O O O' ^ COOCSJ CO in vO CO O ^ h“ 
rocNj *— -i •— 1 12 

ocoocooooooocooocoooocooooooooooo 



OOOOOOOOOOOOOO 00 00 00 0000 0000 00 oooooooooooooo 000 



'^^Ornrnfnfnrnrnrri^ 

oooooooooooooo 

1 

OOCOQOOOOOCQOO 

>i*rnoinr^cocNJc\i^rooocovrc> 

rONOr^mmmf^oO'^’rovffnsO 

omlnln^Oco<^i^Olnmoc^Jrornl— 

(Nj o c\j r -4 CO vf CO m (Nj o o Nf > 
fo ^ o o f\ih- (M o sO c\j 

r-^<Ninini^in^c\jrorornco^r\i ooooooooooooooooooooooooooooooooo 



oooooooooooooo coooocoooooooooooooooo 00 00 00 000 00 

I ) I 



(r^(^rnrnrnrn(^(r^(r^(r^ro<^(^(^ rOforn(r^<r^rO(ri(r^fn(r^fr^(r^(r^rn(^(^fnrnrn(^(r^(r^(r^(r^(rj(r^(r^(r^(r^r^ 

oooooooooooooo oocooo 000 0000 000 00 00 00 000 c 0000000 



OOOQOOOOCOOOOQ OOOOOQCiOOOQOCCiCOOQQQCOOCCC 
h-vfco^Jcor\jO(NJ^Oo— ^' 4 “CM »— •{NjNOOvjf^i'nominLnsOin^'OOsjm 

co<^0'^ovjog>t000^co >oc\j mmLraruninininmc\iroOsOcoocD(^iriror>jr^ooOOr^sO 
fTicooo^in^c\jsoininr^CMnKcccooococooocococcoococccoa'Oc>ooococooooo 
•-^ONtoOsoincoh-'.'MC^oojOinxooa'OOoooo'ooa'Ooooooooa'Ooooo 
Oinrn.— ^OnJ(^Jrns^•ln^^•^s^•(^J<^J r-l,— 11 — 41 — 

c\j vjj* 00 ro ro co ro ro ro ro ro ro •—4 *—4 r— t #— h ^ 1—4 »— 4 1—4 (—4 1—4 *—4 (—4 *—4 »— < 1—4 *—4 ,_4 ,_i *— • f— 4 f »4 



CQOQOCCi 
^inc 'O ococc 
>Orvj •— un o 0 (Nj 
OCOOOOOO 

ooocx^oc 



4 (\J^ 



-•(NJCNJ 



oooooooooooooo OOOCOOOOOOOOOOOCOOOCOOOOOOOOOCOOO 



(NjrvJCM*— 4 C\J (\J (\J •— 4 f— 4 I— 4 ^ 

0000 COOOOOOO 



0000 

men O' > 4 “ 

*— 4 ' 4 ' C\IO 
h-movjm 
•— 4 1—4 r*» 



qooqooqo ^ 

— 40 m o 4 — 

in o ro •— < o ' 4 " Nt w- 

^OOCOO^Or-) —4 

cL 

(\j CO in sT <r sO 12 
o *— » •— 4 ^ o h- in ro r -4 o o 



000 00 000 000 000 00 00 00 00000 00c 0000 



OOOOOOOOOOOOOO o 0000 00 00 00 oooooooooooooo 0000 000c 



>3“in CO o o — 4 (\im '4’in sor^o— ^njrr', N^insor^cocT'O^ 

OOC'O'OCT'OOOCOOOOO 



cNjcOvtinsOh-ccoo— 4 (\jrO' 4 'insOf^cooo— ^c\irn 
^ ^ ^ ^ ^ ^ ^ ^ r\i 04 (XI (\j (\J (Nj r\J c\j (XI c\j ro rn rn ro 



114 



rorOf<^rororof^rorO(Nrnrornrornrororornrornro(N 

OOOOOOOOOOOOOOOOOOOOOOO 

QQQQQOOQCOQQQaQQCQQQQOO 

f^i^oosCo^o^<NooG^>^'rovDc\imr-roaDoooo<fmCMr\ 

<Nj O' c^ o o CO o O ' 00 vO ^ ^ LA >4“ o ro 

o so AO •-^A.r-^o<fNr'^oc^JcoA^^*oo'^Jooo' 
c\j f\j (N CO CO O O CM rn (M ^ rvj •— t O A CM cNj nJ- CO o <f 

oo'^A^Ar*^r^co^rM>4-LnA«— to^^o^o<^JO'0'^<^ 

ooooooooo»-Hr-4rnro<f>f>j'v4-NfoorOvrNj-m*-»f\irnrornrorM^O'0000000000000000 



OOOOCOOCOOOOOOOOOOOOOOOOOOOOCCOOOOOOCOOOOOOOOOOC 



ror<^rnroroc<",rofrrororn<^troro(Nrn 

oooooooooooooooo 

QQCOQOOQQOQOOQOQ 
(NI(NOOOO-4*COOv^f'^0'0 0'C'0'f'^ 
%0 A>tCMcX)^>d“^OrM^fM 

AOOO'COCTO'OAh-Nd-mrMACOf^ 
0>J"^o^'OC'rocoO'0'^ro>4-ooco 
h-OAAoovOco-- tA^vrr^r^a'ogo 

O o o O o O O O O o O O o O O o O O O O O O O o C O O O o O O O A Cvj ^ ^ (N) ^ ^ (N CvJ ro CO m vj- A O' 



OOOOOOOOOOOOOOCOOOOOOOOOO OOOOOOOOOOOOOOOOOOOOOOO 



rorn(^nrn{^rnrornrnroc^rororornrornrnct^ft^(r)rnrnrnrororOrornrnrorr)rornrororororornrnroronrn 
oooooooooooooooo OOOOOOOOOOOOOOOO oooooooooooooooo 

QQOOQOOQOOQOOOCOQQOOOOQQCjOCQOOQOqOCjQQCQQQQCOCCOQ 

ro CO A o m m A O' nO o I— t o 'O >4" O' cvj o CM vO o o \0 ^ o CM ^ ^ CO (N m sC A >4" <f CO CO >d" 'O o nO 

Ocovoor^o^ro vor^cM A^4-C^0' oogcMco^ Acoc\jcoco(\if— ii^ror^'4T\)Oo' CONOCO >4“>4-C'>4-r-^<fNCrococc 
OOcor^AC'roA<fvOAroA(Njr^o^cosOcooOQO^O'OAcoAvOCMCMcoroocovOh-0"Oco'«ornr''O^^C'0 

OOOOOOOOO0'^-H<f^0O'--'^n^4“AOf<^0'CA^^0'•—^^4'^0C^---'^c^^AG'^^0^*<^0'O*--t(^JvfA^^O'4■O1—^ 
^ ^ ^ ^ ^ ^ (\J (\J (\j O ^ ^ ^ ^ C\J (NJ (\| (\| O ^ 04 o ^ ^ ^ ^ (N c\j CM (N ^ CM <N (\I ro *--t (N rg (N og (\| 04 CO ro ^ (\J 



OO OOOOOOOOOOOOOOOOOOOOOOO OOOOOOOOOOOOOOOOOOOOOOO 



sC A A A A>4“ ^fA'O AACO>4"g3 A A AA AA AA43 A A A AA AA A AA A AAA A 
oooooooooooooooo OOOOOOOOOOOOOOOOOOOOOOO 

QOaQOaQQOQQQQQOOQQQQOOOQQOQOQQQQOCQQQOQ 
O' ^ O' CO s4- >4- CO CO o M- o A CO O' o <f ^ CO CO eg ^ o m 43 A o O' O' vJ* 1-4 CO A 
O' CO A O ^ CO CM O 4D o O CO CM O nT O' vC CM O r- 00 sO 43 CO eg ^ A ro A r- 4D 43 CO >4“ v4" CM CM CO CO 
43 r*^ A O ^ 43 O O CO A A <f ^ (M ^ rg CO O ro A A CO ^ CO CM O ^ CO f'- CO ^ 43 CO A CO o O 
O' r- 43 GO >4“ o^ O' CO r- o ro ^ A O" (M ^ CO n 4- sj- »-H A o CM ^ CM n 4- 43 m CO r- CM (M CO r-^ r- r- <f 
04 N4“f'^43 43r^Acn>4TOArOOr-fMO'43^43N^04 0' A0go43CM A43 430 43 M->4-43 43Nf O'O 
OOOOOOOOOr-^O'ArO'— •>4'M“O'^^Ar^43^^C0A<fO404^CMrn^A04r4rOrOrnv4"'4'AAA43r^^CM 



OOOOOOOOOOOOOOOOOOOOOOOOOOOOOCO<^OC OOOOOOOOOOOOCO 

III I 



<fA43r^COCT‘O^CMrnM“A43r^OOO'Of-404rn<fA43^^COO'0^04Cn>4-AsOr^COO'0^04rO<fA4?O^COO'0^ 
m CO ro m m <f <f M" ■ 4" M" M" M" 4" A A A A A A A A A A 43 43 43 43 43 43 43 43 43 43 CO 00 



115 



OOOOOCOOOOOOCOOOOOOOOOOCOO 

ooooooooooocc oooc ooooooooo 



rnroro 

ooo 

OOQ 

COh-CM 

ootnO 
rnCM— < 

r-^^(MOOOOOOOOOOOOOOOOOOOOOOO 



OOOOOOOOOOOOOOOOOOOOOOOOOO 



r 0 r 0 r 0 r 0 (^r 0 r 0 (^<^fnt*nr 0 rt^rnr 0 rnroc 0 rrirnr 0 r 0 <^,rnrn ro 

OOOOOOOOCOCOOOOOOCOOOOOOCO 



OOQCCqQC,QQC.QCjQQC:QQoQQOQOQO 
rn (Nj r- cn r-* ro ^ in nO (N o CO CO ^ (\i ro <\i c\i 

h-vO^COCONO(\JC\ir^rO<^fOh-f\l^rO'^’inr^O^CO< 7 '' 4 “ 0 ^r— 
sOro^ 0 ^f^O(\icoO>tro-^c\i-«^O^^^corO’-^N^’’-^in— ^com 
o (\| ro ^0 O >j- O (M in h~ <M O ro LO vC CO r-» 
nj nj ro (\J nj rvj (\j <N r\j ro ro eg rvj c j f\j ro (\i c\i (\i r\i r\i (\j (\i c\i 



OOOOCOOOOOOOOOOOOOOOOOOOOO 



irur.in 

ooo 



Ooo 

rO^*“ nj 
OOON 
rOin>t 
> 4 -U^nO 
^sonj 

(MCMroOOOOOOOOOOOOOOOOOOOOOOO 

OOOOOOOOOOOOOOOOOOOOOOOOOO 



rMrovfinor^GCC^Of— inirt-^vTin 

oocoaocococococoo'o^crooa^ 



vor^oocT'O^njrOvrmsOf^ 

oooa*oooooooo 



116 



APPENDIX D 



CALCULATION OF ROTOR ELEMENT FLOW ANGLES 



The following is a brief synopsis of the procedure 
contained in Ref. 13 for calculating the outlet relative 
flow angles in a rotor element from the given inlet relative 
flow angle and blade solidity. The reader is referred to 
Ref. 13, Chapter VI, for specific details of low speed 
correlation data. 



As stated in Section III. A, uniform flow conditions at 
the rotor blade edges were assumed. This assumption coupled 
with knowledge of the mass flow rate and rotational speed, 
enables one to calculate the inlet relative flow angle, 
as shown in Fig 18. 

From blade geometry information, the blade solidity, 0" , 



r= ^ 



( 1 ) 



is obtained. At this point, p, , and are given and one may 
calculate , the rotor outlet relative flow angle from 
correlation curves depicted in Ref 13. The equation used to 
determine jS^ , is the following. 






( 2 ) 



where K is the angle between the tangent to tne blade mean 
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camber line and the axial direction (Fig 18) . This is 
obtained from the blade geometry data. ^ is the low speed 
deviation angle which is obtained from the correlation 
curves in Ref. 13. The following equations show the 

relationship between S and the correlation data. 

1= ,3) 










The variables m,K$) ,K$), and ) , are all values 

5K t <«> 

waich are obtained from the correlation curves and are all 
functions of the given blade geometry. The quantity, <|> , is 
the blade camber angle and again is obtained from the blade 
geometry data. Once all the variables are obtained from the 
correlation data, equation (4) is solved for the deviation 
angle for an uncambered blade section, and then equation 
(3) is solved for the deviation angle, ^ . One now 

calculates ^ 2 . ftom equation (2) for the blade element. With 
^ 3 . now a known quantity, one now calculates the absolute 
flow angle, c(^ , from uniform flow assum-ptions. 



An example follows for no-de numbers 43 and 57 (Fig 5). 
From Ref. [12], Table II, the following quantities are 
obtained assuming the angle of incidence, i, (Fig 13) is 
zero and therefore the inlet relative flow angle, p, , is 
equal to . 

p, = k, = 61.38 
r = 1.31162 






6 . 95 






1 ta 



=3 



= 54.93* 

jt \ 

~ET ) tnax = 0.035 
tip radius = 18.25 in 
hub radius = 9.125 in 



Assniraing" uniform flow at the rotor inlet and a 
rotational speed of" 4359.5 HPM, the following quantities are 
determined from the rotor inlet velocity diagram (Fig IB). 

(107. (r iWwi/sec) (l4-4- iwVFt*) 

V Qi = ^ A - ^ ; 

T (o.Oft TT (ie.2S*- i.125.') IN* 

Vm = 246.802 Ft/sec 

where the area A, is determined from the hub and tip radii 
and the density is assumed to be 0.08 Ibm/cu ft. 



Now one is ready to obtain the correlation data. From 
Ref. [13], Fig 162, with 3, = 61. 88*^ and T = 1.3062, 



2.50* 

Ki © 

From Ref. [13], Fig 162, with = 61.38 and (T =1.3052 , 
m = 0.235 

From Ref. [13], Fig 172, with t/c) max = 0.0350, 

t = 0.29 

From Ref. [13], page 222, one uses the following value of 
(K^) sh for 65-series blades. 



kS) sh = 1.0 

At this point all the necessary data has been obtained for 
equations (3) and (4) , 

C = 29) (2. 50) = 0.725* 

From equation (3) , 

S = 0.725*f 0.235 (6.95'’) = 2.36* 

Finally, equation (2) gives the desired value of , 



1 19 











= 54.93'’ 


+ 2 


.36*= 57 


. 29* 






At 


this 


point 


the 


relati ve 


flow angle 


for node 57 


has been 


obta 


ined , 


57. 


29^ . 


These 


two 


values 


of 


relat 


ive flow 


angl 


es. 


ji, = 61 


. 88 * 


for node 


: 43 


and = 


57. 29 


for 


node 57 , 


are 


then 


read in 


the 


program 


as 


input 


data 


for 


numerical 


comp 


utat ion. 


















This 


process 


is 


repeated 


at 


each req 


uired 


blade 


element 



section for the proper outlet relative flow angle. 
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Figure 18 - NOMENCLATURE FOR CASCADE BLADE 
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